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Abstract 

When a black hole evaporates, there arises a net energy flow from the black hole into its outside 
environment due to the Hawking radiation and the energy accretion onto black hole. Exactly speaking, 
any thermal equilibrium state has no energy flow, and therefore the black hole evaporation is a 
nonequilibrium process. To study details of evaporation process, nonequilibrium effects of the net 
energy flow should be taken into account. The nonequilibrium nature of black hole evaporation is 
a challenging topic which includes not only black hole physics but also nonequilibrium physics. In 
this article we simplify the situation so that the Hawking radiation consists of non-self-interacting 
massless matter fields and also the energy accretion onto the black hole consists of the same fields. 
Then we find that the nonequilibrium nature of black hole evaporation is described by a nonequilibrium 
state of that field, and we formulate nonequilibrium thermodynamics of non-self-intcracting massless 
fields. By applying it to black hole evaporation, foUowings are shown: (1) Nonequilibrium effects of the 
energy flow tends to accelerate the black hole evaporation, and, consequently, a specific nonequilibrium 
phenomenon of semi-classical black hole evaporation is suggested. Furthermore a suggestion about 
the end state of quantum size black hole evaporation is proposed in the context of information loss 
paradox. (2) Negative heat capacity of black hole is the physical essence of the generalized second 
law of black hole thermodynamics, and self-entropy production inside the matter around black hole is 
not necessary to ensure the generalized second law. Furthermore a lower bound for total entropy at 
the end of black hole evaporation is given. A relation of the lower bound with the so-called covariant 
entropy bound conjecture is interesting but left as an open issue. 
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1 Introduction 

Black hole evaporation is one of interesting phenomena in black hole physics [T] . A direct treatment of time evolution 
of the evaporation process suffers from mathematical and conceptual difficulties; the mathematical one will be seen 
in the dynamical Einstein equation in which the source of gravity may be a quantum expectation value of stress- 
energy tensor of Hawking radiation, and the conceptual one will be seen in the definition of dynamical black hole 
horizon. Therefore an approach based on the black hole thermodynamics is useful [2] [^ . 

Exactly speaking, dynamical evolution of any system is a nonequilibrium process. If and only if thermodynamic 
state of the system under consideration passes near equilibrium states during its evolution, its dynamics can be 
treated by an approximate method, the so-called quasi-static process. In this approximation, it is assumed that 
the thermodynamic state of the system evolves on a path lying in the state space which consists of only thermal 
equilibrium states, and the time evolution is described by a succession of different equilibrium states. However once 
the system comes far from equilibrium, the quasi-static approximation breaks down. In that case a nonequilibrium 
thermodynamic approach is necessary. For dissipative systems, the heat flow inside the system can quantify the 
degree of nonequilibrium nature [H [51 [SI [I] ■ 

For the black hole evaporation, when its horizon scale is larger than Planck size, it is relevant to describe the black 
hole itself by equilibrium solutions of Einstein equation, Schwarzschild, Reissner-Nortstrom and Kerr black holes, 
because the evaporation proceeds extremely slowly and those equilibrium solutions are stable under gravitational 
perturbations [8]. The slow evolution is understandable by the Hawking temperature [T] which is regarded as an 
equilibrium temperature of black hole, 



where M is the black hole mass, rripi is the Planck mass and the units c — h = ks = 1 and G = l/rripi are used. 
Obviously a classical size black hole (M ^ mpi) has a very low temperature. This means a very weak energy emission 
rate by the Hawking radiation which is proportional to (2GM)^ Tg due to the Stefan-Boltzmann law. Therefore the 
quasi-static approximation works well for the black hole itself during its evaporation process. However the outside 
environment around black hole may not be described by the quasi-static approximation because of the energy flow 
due to the Hawking radiation. The Hawking radiation causes an energy flow in the outside environment, and that 
energy flow drives the outside environment out of equilibrium. As indicated by Eq. (jl.ip , the black hole temperature 
and the energy emission rate by black hole increase as M decreases along the evaporation. The stronger the energy 
emission, the more distant from equilibrium the outside environment. Therefore the nonequilibrium nature of 
the outside environment becomes stronger as the black hole evaporation proceeds. At the same time, the quasi- 
static approximation is applicable to the black hole itself since equilibrium black hole solutions are stable under 
gravitational perturbation. Hence, in studying detail of evaporation process, while the black hole itself is described 
by quasi-static approximation, but the nonequilibrium effects of the energy flow in the outside environment should 
be taken into account. 

In the above paragraph, the energy accretion onto black hole is ignored. However if the temperature of outside 
environment is non-zero and lower enough than the black hole temperature, then the black hole evaporates under the 
effect of energy exchange due to the Hawking radiation and the energy accretion. In this case the same consideration 
explained above holds and we recognize the importance of the net energy flow from black hole to outside environment. 
Dynamical behaviors of black hole evaporation will be described well by taking nonequilibrium nature of the net 
energy flow into account. 



2 



In SecIJl we introduce a simple model of black hole evaporation to examine the net energy flow in the outside en- 
vironment, where the matter fields of Hawking radiation and energy accretion are represented by non-self- interacting 
massless fields for simplicity. Secl3]is devoted to construction of nonequilibrium thermodynamics of that field. Then 
it is applied to the black hole evaporation. Secl4]reveals that the nonequilibrium effect tends to accelerate the evap- 
oration process and, consequently, gives a suggestion about the end state of quantum size black hole evaporation 
in the context of the information loss paradox. Sec\E\ reveals that the generalized second law is guaranteed not by 
self-interactions of matter fields around black hole which cause self-production of entropy inside the matters, but 
by the self-gravitational effect of black hole appearing as its negative heat capacity in Eq. (|2.2p . Readers can read 
SecsU] and [5] separately, and may skip over Sec|4]to see discussions on generalized second law in SeclS] Finally Secl6] 
concludes this article with comments for future direction of this study. 

Throughout this article except for Eq. ()l.ip . Planck units c — h — G — ks — I are used. Then the Stefan- 
Boltzmann constant becomes a := n^/GO which is appropriate for photon gas. When one consider non-self- 
interacting massless matter fields, as indicated in Secl31 it is necessary to replace cr by its generalization, 

Ntt^ N , , 

where N :— nt, + (7/8) n/. Here rtf, is the number of inner states of massless bosonic fields and Uf is that of massless 
fermionic fields, {rih — 2 for photons.) Furthermore, at least when the black hole temperature is lower than 1 
TeV (upper limit by present accelerator experiments), it is appropriate to estimate the order of N by the standard 
particles (inner states of quarks, leptons and gauge particles of four fundamental interactions), 

TV -100. (1.3) 

This denotes a' ~ 10. Throughout this article, we simply assume that N is independent of black hole temperature 
and this estimate (|1.3p holds always for semi-classical black hole evaporation {Tg < 10^^ TeV). 



2 Thermodynamic model of black hole evaporation 

According to the black hole thermodynamics [Dill [3], a stationary black hole is regarded as an object in thermal 
equilibrium, a black body. For simplicity, let us consider a Schwarzschild black hole. Its equations of states as a 
black body are 

^o = ^— = — , S.^^-^^ttR^, (2.1) 
^ 8TTTg 2 ' ^ 167rr2 f ^ ' 

where Eg, Rg, Tg, Sg correspond respectively to mass energy, horizon radius. Hawking temperature and Bekenstein- 
Hawking entropy. Obviously the radius Rg decreases when this body loses its energy Eg. The black hole evaporation 
is represented by the energy loss of this black body. 
The heat capacity Cg of this body is negative. 

The negative heat capacity is a peculiar property of self-gravitating systems [Q, . Therefore the energy Eg includes 
self-gravitational effects of a black hole on its own thermodynamic state. Furthermore it has already been revealed 
that, using the Euclidean path-integral method for a black hole spacetime and matter fields on it, an equilibrium 
entropy of whole gravitational field on a black hole spacetime is given by the Bekenstein-Hawking entropy [lOj . 
This means Sg in Ea. (|2.ip is the equilibrium entropy of whole gravitational field on black hole spacetime, and the 
gravitational entropy vanishes if there is no black hole horizon. Hence we find that energetic and entropic properties 
of a black hole are encoded in the equations of states ()2.ip . Hereafter we call this black body the black hole. 

As mentioned SeclU the nonequilibrium nature of black hole evaporation arises in the matter fields around 
black hole due to the net energy flow by Hawking radiation and energy accretion onto the black hole. When we 
consider arbitrary dissipative matter fields as the Hawking radiation and energy accretion, we immediately face a 
very difhcult problem how to construct a nonequilibrium thermodynamics for arbitrary dissipative matters. This is 
one of the most difficult subjects in physics [11 [5l O [7] . To avoid such a difhcult problem and for simplicity, let us 
consider non-sclf-intcracting massless fields to represent the Hawking radiation and energy accretion. For example, 
photon, graviton, neutrino (if it is massless) and free Klein-Gordon field (□$ = 0) are candidates of such matter 
fields, and they possess the generalized Stefan-Boltzmann constant a' given in Eq. (|1.2p . Hereafter we call these 
fields the radiation fields. 

As mentioned above, nonequilibrium phenomenon is one of the most difficult subjects in physics. It is impossible 
at present to treat the nonequilibrium nature of black hole evaporation in a full general relativistic framework. Hence 
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Figure 1: NE model. Left panel shows energy flow between black hole and heat bath. Right one shows some variables 
and particle orbits of radiation fields. The radiation fields sandwiched by black bodies are in two-temperature steady 
state. 

we resort to a simplified model to examine the nonequilibrium effects of net energy flow in the outside environment 
around black hole [TT]: 

Nonequilibrium Evaporation (NE) model: Put a spherical black body of temperature Tg in a heat bath of 
temperature Th{< Tg), where the equations of states of the spherical black body is given by Eq. (|2.ip and 
we call the black body the black hole. Let the heat bath (the outer black body of temperature T/i) be made 
of ordinary materials of positive heat capacity. Then hollow a spherical region out of the heat bath around 
black hole as seen in Fig[TJ The hollow region is a shell-like region which is concentric with the black hole 
and separates the black hole from the heat bath. This region is filled with matter fields emitted by black hole 
and heat bath. Those matter fields are the non-self-interacting massless fields possessing generalized Stefan- 
Boltzmann constant a' given in Eq. (|1.2p . and we call these fields the radiation fields. This model consists 
of three parts, black hole, heat bath and radiation fields. Furthermore we consider the case that the whole 
system is isolated and the total energy of the three parts is conserved. 

The isolated condition reflects the whole universe including an evaporating black hole. The temperature difference 
(Tg > Th) causes a net energy flow from the black hole to the heat bath. This energy flow drives a relaxation process 
of the whole system because of the isolated condition. For ordinary systems of positive heat capacity, relaxation 
process reaches an equilibrium state at the end of the process. But the relaxation in NE model may not reach a 
total equilibrium state of the whole system, because the temperature difference Tg — Th may increase along the 
decrease of Eg due to the negative heat capacity Cg given in Eq. (|2.2p . This causes the decrease of Rg due to the 
equations of states p.ip . Therefore the relaxation process in NE model corresponds to the black hole evaporation. 

Here let us comment about terminology. There may be an objection that the term "relaxation" is not suitable 
for the case of increasing temperature difference. But in this article, please understand it means the time evolution 
arising in isolated inhomogeneous systems. 

If a full general relativistic treatment is possible, the Hawking radiation experiences the curvature scattering to 
form a spacetime region filled with interacting matters and some fraction of Hawking radiation is radiated back to 
the black hole from that region. The heat bath in the NE model is understood as a simple representation of not 
only matters like accretion disk but also such region formed by curvature scattering. 

Here, as an objection to the NE model, one may remember the so-called Tolman factor which appears in the 
"equilibrium" temperature of radiation fields around black hole: Exactly speaking, equilibrium of any matter fields 
around black hole is "local" equilibrium. It has already been known that, when the radiation fields around black 
hole are in local equilibrium with the black hole, the local equilibrium temperature Teq{r) of radiation fields (not 
of black hole) at a spacetime point of areal radius r from the center of black hole is T^g[r) = Tgj ^\ — Rg/r, where 
Tg is the Hawking temperature and Rg is the horizon radius. This Teq{r) is obtained in a full general relativistic 
framework of equilibrimn thermodynamics, and the factor — Rg/r is the Tolman factor T2j. One may think 

that, because T^q —> cxd as r i?g, it is unreasonable to assign Tg to black hole as its temperature. But let us 
emphasize that Tf,q{r) is not black hole temperature but the equilibrium temperature of radiation fields. We can 
understand Teq{r) as follows: In order to retain equilibrium of radiation fields against external gravitational force 
by black hole, a temperature of radiation fields higher than the asymptotic value Tg is required, since the higher 
temperature denotes the higher pressure against external gravitational force. The Tolman factor describes the effect 
of external gravity on the equilibrium radiation fields, and becomes unity if the external gravity vanishes. The local 
equilibrium temperature of radiation fields may count an "intrinsic" temperature of black hole and an additional 
gravitational effect in Tolman factor. It may be reasonable to regard the asymptotic value Tg as an intrinsic black 
hole temperature. Hence we assign Tg to the black hole in NE model. But it is turc that the NE model is not 
a full general relativistic model, and ignoring, for example, gravitational rcdshift and curvature scattering on the 
"nonequilibrium" radiation fields propagating in the hollow region. Although the NE model may be too simple, let 
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us try to investigate nonequilibrium effects of the net energy flow from black hole to its outside environment in the 
framework of NE model. 

Furthermore, we put the quasi-equilibrium assumption to utilize the "equilibrium" equations of states (j2.ip for 
black hole, and the fast propagation assumption to treat the nonequilibrium nature as simple as possible: 

Quasi-equilibrium assumption: Time evolution in the NE model is not so fast that the evolutions of black 
hole and heat bath are approximated well by the quasi-static process individually, while the quasi-static 
approximation is not valid for the radiation fields due to the net energy flow from black hole to heat bath. 
Then, it is reasonable to use Ea. (|2.ip as the equations of states for black hole. Furthermore, since Schwarzschild 
black hole is not a quantum one, following relation is required, 

i?9 > 1 • (2.3) 

Fast propagation assumption: The volume of hollow region is not so large that particles of radiation fields travel 
very quickly across the hollow region. Then the retarded effect on radiation fields during propagating in the 
hollow region is ignored. 

There are two points which we should note here. The first point is about the thermodynamic states of black hole 
and heat bath under the quasi-equilibrium assumption. This denotes the temperatures Tg and Th are given by 
equilibrium temperatures at each moment of evaporation process. Therefore we can regard Tg and Th as constants 
within a time scale that one particle of radiation fields travels in the hollow region until absorbed by black hole or 
heat bath. This is consistent with the fast propagation assumption. 

The second point is about the thermodynamic state of radiation fields. In the hollow region, the radiation fields 
of different temperatures Tg and T^ are simply superposed, since the radiation fields are of non-self-interacting 
(coUisionless particles gas). This means the radiation fields are in a two-temperature nonequilibrium state. Fur- 
thermore, because Tg and Th are constant while one particle of radiation fields travel across the hollow region, it is 
reasonable to consider that the radiation fields have a stationary energy flow from black hole to heat bath within 
that time scale. Hence, at each moment of time evolution of NE model, the thermodynamic state of radiation fields 
is well approximated to a macroscopically stationary nonequilibrium state, which we call the steady state hereafter. 
Consequently, time evolution of radiation fields is described by the quasi-steady process in which the thermodynamic 
state of radiation fields evolves on a path lying in the state space which consists of steady states, and the time 
evolution is described by a succession of different steady states. Therefore we need a thermodynamic formalism of 
two-temperature steady states for radiation fields. The steady state thermodynamics for radiation fields has already 
been formulated in |13j . which is summarized in next section. 

3 Steady state thermodynamics for radiation fields 

Before proceeding to the black hole evaporation, two-temperature steady state thermodynamics for radiation 
fields [13] is explained in this section. Subsec lSTTl introduces the minimum tools required to apply to the NE 
model. Readers may skip over Subsec i3?2l to read remaining Sees HI [5] and [H Furthermore, since SecsH] and[5] are 
written separately, one can also skip over SecHlto see the generalized second law in SeclH 

Subsecl321 exhibits a more detail of two-temperature steady state thermodynamics for radiation fields. Although 
a full understanding of it is not necessary for black hole evaporation, but it may be helpful to understand, for 
example, the free streaming in the universe like cosmic microwave background and/or the radiative energy transfer 
inside a star and among stellar objects. Keeping future expectation of such applications in mind, Subsec l5?^ is 
placed here. 

3.1 Minimum tools for applying to NE model 

To concentrate on investigating two-temperature steady states of radiation fields, we consider a model which can 
be realized in laboratory experiments. According to the NE model, we introduce the following model named SST 
after the Steady State Thermodynamics. Then, after constructing the steady state thermodynamics for radiation 
fields, we will modify the SST model to the NE model in SecsH] and O 

SST model: Make a vacuum cavity in a large black body of temperature T^ut and put an another smaller black 
body of temperature Tin (t^ Tout) in the cavity as seen in Figj^] For the case T^n > T^^t, the radiation fields 
emitted by two black bodies causes a net energy flow from the inner black body to the outer one. When 
the outer black body is isolated from outside world and the heat capacities of two black bodies are positive 
definite, the whole system which consists of two black bodies and radiation fields between them relaxes to a 
total equilibrium state in which two black bodies and radiation fields have the same equilibrium temperature. 
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Figure 2: SST model. Radiation fields between outer and inner black bodies are in steady state. 



It should be emphasized that Tin and Tout are "equilibrium" temperatures of black bodies, respectively. Outer 
body is in an equilibrium state, and inner one is also. But their equilibrium sates are different from each other 
since Ti„ ^ Tout- Then the radiation fields between them is in a nonequilibrium state. By keeping temperatures Tin 
and Tout constant, the nonequilibrium state of radiation fields becomes a macroscopically stationary nonequilibrium 
state, which we call a steady state. 

The difference of SST model from the NE model is that heat capacities of two black bodies are positive definite 
and the system does not necessarily has symmetric geometry. Because of the positive heat capacity, the relaxation 
process in the SST model leads the whole system to a total equilibrium state, while the black hole evaporation 
increases the temperature difference between black hole and heat bath due to the negative heat capacity of black 
hole. But as for the NE model, the radiation fields in the SST model are also non-self-interacting massless fields. 
The SST model ignores gravitational interactions among two bodies and radiation fields. Furthermore we consider 
the case satisfying following two conditions according to the quasi-equilibrium and fast propagation assumptions; 
the first is that the evolution of each black body is of quasi-static during the relaxation process of the whole system, 
and the second condition is that the volume of cavity is so small that retarded effect on radiation fields is ignored. 
Then, as discussed SeclSJ the evolution of radiation fields is of quasi-steady. Due to the Stefan-Boltzmann law, each 
steady state composing the quasi-steady process of the radiation fields has a net energy fiow 

J sst = O- {T^n - Tout ) ^in , (3-1) 

where Ain is the surface area of inner black body. This Jsst equals the net energy exchanged par a unit time between 
the two black bodies via the radiation fields. 

A consistent thermodynamic framework for the steady states of radiation fields has already been constructed [13j . 
The outline of the construction of steady state internal energy and steady state entropy is as follows. 

For the first we consider a massless bosonic gas. The energy of the gas is given by 

/dp^ 3 

where p is a momentum of particle, a; is a spatial point, e^g is an energy of particle of p at x, and gp^g and dj^^g are 
respectively the number of states and the average number of particles at a point (p, x) in the phase space of the 
gas. For the entropy, we refer to the Landau-Lifshitz type definition of nonequilibrium entropy Sb for bosonic gas 
(see §55 in chapter 5 of [14]), 

f dp^ 

^6 := y (2 7r)3 ^^'^ ^^'^ [ ( 1 + ^P^s ) In ( 1 + dft^g ) - dp^g In dp^g ] . (3.3) 

It has also been shown in §55 of [M] that the maximization of Sb for an isolated system {SSb = 0) gives the 
equilibrium Bose distribution. This is frequently referred in many works on nonequilibrium systems as the H- 
theorem. However in §55 of [14] . concrete forms of gj;g and dj^g are not specified, since an arbitrary system is 
considered. 

In the SST model, the radiation fields are sandwiched by two black bodies, and its particle is massless and 
"coUisionless" . Then we can determine Cp'^j, gp^g and dpg to be 

f-p,x = w , gp^x = rib , dp^s = r ,rJf^ 7 > (3-4) 

exp[uj/T{p,x)]~l 

where the frequency of a particle a; = |p| , n;, is the number of inner states of bosonic gas which is assumed to be 
constant as mentioned in Eq. (|1.3p . and T{p, x) is given by 

T{p,x):=r "' ^"^^^ ""^^ , (3.5) 

for p = Pout at X 
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where pi„ is momentum of a particle emitted by the inner black body and Pout emitted by the outer one. The 
x-dependence of T(p, x) arises, because the directions in which particles of Pm and Pout can come to a point x vary 
from point to point. 

From the above, we obtain the steady state energy and entropy, 

Eb = J dx^Cbix) , eb{x) := iab {gin{x)Tif^ + gout{x)T^ut) (3.6a) 

Sb^ J dx^Sbix) , SbiS) := ^ ( 5,„(f) T^l + gout{x) T^lt ) , (3.6b) 

where integrals in Eq. (|3.1ip are used, fjf, :— rif,7r^/120, gin{x) is the solid angle (divided by An) covered by directions 
of Pin at x as shown in FiglJl and gout{x) is similarly defined with Pout- By definition, we have gin{x) + gout{x) = 1. 
On the other hand, if the radiation fields are in equilibrium of temperature T, the ordinary equilibrium thermo- 
dynamics determine equilibrium energy density e and entropy density s to be e = 4cr;,T^ and s = (16crh/3)T^. 
We find that the steady state energy and entropy are given by a simple linear combination of the values which 
are calculated as if the radiation fields are in equilibrium of temperature Ti„ or Tout- This is consistent with the 
collisionless nature of radiation fields. 

Next we consider a massless fermionic gas, the formula of energy (|3.2p holds for fermions as well. But the 
formula of entropy (|3.3p is replaced by [Mj 

[ dp^ 

Sf := - J ^2Tr)^ '^^^ ^P'^ ii^^ dp.s ) In ( 1 - dj^^s ) + dj^^s In dff.s ] . (3.7) 

Then, following the same procedure as for the massless bosonic gas with replacing the distribution function by 
dp,x = [exp[a;/T(p, x)] + we can obtain the steady state energy and entropy of a massless fermionic gas, 

Ef ^ --^Eb , Sf = --^Sb, (3.8) 

a Ub o Ub 

where integrals in Eq. (|3.1ip are used, and ?t./ is the number of inner states of fermionic gas which is assumed to 
be constant as mentioned in Ea. (|1.3p . Hence, when the black bodies in SST model emit n/ massless fermionic 
modes and massless bosonic modes, the steady state energy Eg^t and entropy Sgst of radiation fields are given 
by Eas. (|3.6p with replacing ab by cr' given in Eq.pT 



Esst ^ I dx^Csstix) , Csstix) -.^ Aa' {gi„{x)Ti'^^+ goutix)T^ut) {3.9a) 

16 a' 



Ssst^ J dx^Ssstix) , Ssst(^) := — ^ (g^n{x)T^l + gout{x)T^^t) , (3.9b) 

At equilibrium limit Tin — Tout, these densities Csst and Ssst become equilibrium ones e — Aa'T'^ and s = (16a/3)T'^ 
respectively, where an identical equation gin (a?) + gout {x) = 1 is used. 

By defining the other variables like free energy with somewhat careful discussions, it has already been checked 
that the 0th, 1st, 2nd and 3rd laws of ordinary equilibrium thermodynamics are extended to include two-temperature 
steady states of radiation fields [13]. This means that the steady state thermodynamics for radiation fields has 
already been constructed in a consistent way. Especially on the steady state entropy, the total entropy of the whole 
system which consists of two black bodies and radiation fields increases monotonously during the relaxation process 
of the whole system, 

dS,n + dSout + dSsst>0, (3.10) 

where the equality holds for total equilibrium states, Sin and Sout are entropies of two black bodies which are given 
by ordinary equilibrium thermodynamics, and a simple sum of total entropy Sin + Sout + Ssst is assumed for not 
only equilibrium states but also steady states. This indicates that the second law holds for steady states of radiation 
fields. 



Here for the convenience to follow Eqs. p.6p and (j3.8p . we list useful formulae; 
dx = — / dx- 



t= - 1 15 Jq + 1 120 



x'^ TT^ IItt^ 



dx ln(e^ - 1) = — — / dx ln(e=" + i) = + F 

- 1 ^ ^180 io + 1 ^ ^60 

where F = [(ln2)2 - tt^] (ln2)2/6 -I- 40(4, 1/2) + (71n2/2) ((3). Here C{z) is zeta function and 0(z, s) is modified 
zeta function (Apell's function). 
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3.2 A more detail 



This subsection exhibits a more detail of two-temperature steady state thermodynamics for radiation fields. A 
peculiar point of radiation fields is the collisionless nature of composite particles. Radiation fields are non-dissipative 
matter. A representative of them is a photon gas. Before a construction of steady state thermodynamics for radiation 
fields, there were some existing works on nonequilibrium radiation fields. 

A traditional treatment of radiative energy transfer, for example that in a star [S], has been applied to a mixture 
of radiation fields (photon gas) with a matter like a dense gas or other continuum medium. In such a traditional 
treatment, successive absorption and emission of photons by components of medium matter makes it possible to 
consider radiation fields as if in local equilibrium states whose temperatures equal those of local equilibrium states 
of medium matter. However, when the radiation fields are in "vacuum" space, the idea of local equilibrium is never 
applicable to radiation fields because they are non-dissipative (see §63 in [l4]). 

Here let us look at ordinary dissipative systems very briefly. For ordinary dissipative systems, the so-called 
extended irreversible thermodynamics treats successfully nonequilibrium states [4]. However this is applicable not 
to any highly nonequilibrium state but to a state whose entropy flux is well approximated up to second order in 
the expansion by the heat flux of the nonequilibrium state. On the other hand, a steady state thermodynamics has 
already been suggested for dissipative systems like heat conduction, shear flow, electrical conduction and so on [S]. 
Also a heat flux appears as a consistent state variable in those steady state thermodynamics. However the range 
of its application is limited [51 [7] . Although present nonequilibrium thermodynamics for dissipative systems have 
some restriction on their applicability, the point is that a heat flux plays the role of consistent state variable which 
quantifies a degree of nonequilibrium nature of ordinary dissipative systems. When one deals with nonequilibrium 
dissipative systems, it is usual to place an interest on the heat flux. 

Notable works on nonequilibrium radiation fields in "vacuum" are given by Essex [15]. Here it may be helpful 
to point out that a "heat" arises by dissipation. No heat flux exists in nonequilibrium radiation fields, but an 
energy flux in radiation fields may correspond to a heat flux in dissipative systems. It was natural that Essex 
considered an energy flux in nonequilibrium radiation fields (photon gas) in vacuum. Essex has shown that, contrary 
to the success in extended irreversible thermodynamics, the entropy flux for nonequilibrium radiation fields is 
NOT expressed by the expansion by energy flux of those radiation fields. Even if the same method of extended 
irreversible thermodynamics is applied to nonequilibrium radiation fields, the energy flux becomes inconsistent with 
the nonequilibrium free energy. This inconsistency will be looked over in Eq. (|3.18p later in this section. This means 
the energy flux does not work as a consistent state variable for nonequilibrium radiation fields in vacuum. 

Apart from Essex's works, there were other works on nonequilibrium radiation fields (photon gas) in vacuum 
emitted by nonequilibrium ordinary matters [16j . They used the information theory. The basis of information theory 
is the assumption that nonequilibrium entropy is given by lnd„e, where dm is a nonequilibrium distribution function 
defined case by case according to the system under consideration (see for example |1| in which the information theory 
is also explained) . Applying the information theory to the total system composed by nonequilibrium radiation fields 
and its source matter, complicated distribution functions for general nonequilibrium radiation fields and source 
matter have been suggested in [TOj. However, after those works were reported, it is revealed in [7] that, at least 
for a matter whose components are colliding and interacting with each other, the distribution function for a steady 
state of that matter derived by the information theory does NOT qualitatively agree with that derived by a steady 
state Boltzmann equation. Furthermore it is also concluded in [7] that the nonequilibrium temperature of that 
matter determined by the information theory has no physical meaning. The information theory does not always 
work well. Therefore, because the suggested distribution function of nonequilibrium radiation fields depends on 
nonequilibrium temperature of source matter derived by the information theory, the reliability of the distribution 
function may not be given in 1(5]. There is no confirmed form of distribution function of nonequilibrium radiation 
fields emitted by nonequilibrium matter. Hence, to avoid the difficult problem on nonequilibrium state of source 
matter, we simply assume in the SST model that the source bodies are in equilibrium states. 

From the above, we recognize the following three facts: (1) The traditional treatment of radiative transfer is 
applicable only to a mixture of radiation fields with a matter which is dense enough to ignore the vacuum region 
among components of the matter. (2) Energy flux is not a consistent state variable for nonequilibrium radiation 
flelds in vacuum, and therefore a special approach different from that to dissipative systems is required to understand 
nonequilibrium radiation fields. (3) A consistent thermodynamic formulation for a system including nonequilibrium 
radiation fields in vacuum has not been accomplished, and therefore a consistent nonequilibrium order parameter 
for radiation fields has not been obtained so far. 

At least for two-temperature steady states of radiation fields in vacuum, a thermodynamic formalism is accom- 
plished and a consistent nonequilibrium order parameter for steady states is obtained in [13] in the framework of 
SST model. Exactly speaking, the radiation fields in SST model is in local steady states, because the distribution 
function in Eq. ()3.4|) and its fcrmion version have af-dependence. The radiation fields in a sufficiently small region 
are in a local steady state, but that local steady state may be different from a local steady state in the other small 
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region. Therefore state variables for SST model should be defined as a function of spatial point x. The extensive 
variable is to be understood as a density. 

Let us exhibit consistent steady state variables for radiation fields. See [13] for detail discussions to justify the 
following definitions of state variables. 

Steady state internal energy density esst{x) and entropy density Ssst{x) 

These are already defined in Eq. (|3.9p . 

Steady state Pressure tensor (in 3-dim. space) 

One may naively expect that the pressure of "steady" state is a global quantity, since a pressure gradient in an 
ordinary dissipative system accelerates components of that system to cause a dynamical evolution. However this is 
not true of radiation fields whose particle is "coUisionless" . As seen below, P^f. becomes a function of x because of 
x-dependence of distribution function dp^s- 

In general, pressure is defined by the momentum flux, the amount of momentum carried by composite particles 
par unit area and unit time. For equilibrium states, momentum flux is homogeneous and isotropic, then equilibrium 
pressure becomes a scalar quantity. However for noncquilibrium states, momentum flux is not homogeneous and/or 
isotropic, then the pressure should be defined as a tensor. 



where p and are, respectively, spatial magnitude and components of momentum of a particle of radiation fields. 
Trace of PlL becomes 



At equilibrium limit Tin = Tout, this trace becomes equilibrium pressure, (4CT'/3)r^, where an identical equation 
ginix) + gout{x) = 1 is used. 

Steady state free energy density fsst 

By a requirement that the differential of free energy by volume gives the minus of pressure as for ordinary equilibrium 
thermodynamics, we can obtain 



At equilibrium limit Tin — Tout, this becomes equilibrium free energy, — (4cr'/3)T''*- 
Steady state chemical potential 

Chemical potential in general can be interpreted as a work needed to add one particle to the system under consid- 
eration. Because particles of radiation fields are coUisionless, no work is needed to add a new one into radiation 
fields. This is the case for either equilibrium or steady states. Indeed, the chemical potential of radiation fields 
(photon gas) in equilibrium is zero. Therefore the steady state chemical potential of radiation fields is zero as well. 

Steady state temperature Tsst{x) 

By a requirement that the differential of fg^t by Tgst gives the minus of entropy density {dfsst/dTggt = ~Ssst) as 
for equilibrium ordinary thermodynamics, we can obtain 




(3.12) 




(3.13) 




(3.14) 



(3.15) 



At equilibrium limit Ti 



Tout, this becomes equilibrium temperature. 



Intensive steady state order parameter r 



Energy fiux j{x) is defined by 



j{x) ■■= jn 



(3.16) 
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where n is a unit vector in the direction of total momentum of particles at x, and 



j N J dp^uj dp^s: cos (j) , (3.17) 

where u} = p is the energy (frequency) of a particle of radiation fields, and (j> is the angle between n and p. If j{x) is 
adopted as an intensive steady state order parameter, its conjugate state variable should also be defined well. Such 
a conjugate variable would be defined by the differential of free energy by j{x). However it has shown in [T3] such 
a conjugate variable vanishes, 

^ = 0. (3.18) 
dj 

This means j{x) is not a consistent state variable, since its thermodynamic conjugate variable does not exist. 
Hence, instead of energy flux, we adopt the temperature difference as an intensive steady state order parameter, 

T := r„ - Tout . (3.19) 

This is obviously intensive variable and satisfies a natural requirement r = at equilibrium limit Tin = Tout- This 
T is consistent with the first law and the concavity of free energy as looked over in Eas. (|3.22p and (|3.23p . 

Extensive steady state order parameter density iIj{x) 

After ordinary equilibrium thermodynamics, we define an extensive steady state order parameter as a thermody- 
namic conjugate variable to r using the differential of free energy density. Hence we define as 

^(f) := = 1^ g^„(f) (Tl - Ti,) . (3.20) 

This satisfies a natural requirement ip = ai equilibrium limit Tin = Tout- 
It may be useful to rewrite ip{x) as 

1p{x) = gtn{x) gout{x) (SeqiTin) - Seq{Tout)) , (3-21) 

where Seq(T) = (16CT'/3)r^ is the equilibrium entropy density of radiation fields of temperature T. It seems very 
natural and reasonable that a difference of entropies quantifies a degree of nonequilibrium nature of steady states. 

Prom zeroth to third laws and concavity of free energy 

Zeroth law, the existence of steady states, is the existence of systems which realize steady state radiation fields as 
shown in Fig [21 

First law can be checked from the above definitions of state variables. We can obtain the following equation, 

desstix) ^ ^ = Tsstix) dssst{x) + T dipix) ^ , (3.22) 



af— fixed 



-fixed 



where we fixed x since a local steady state is considered, and consequently a "work term" including pressure does 
not explicitly appear in this relation. A work term will appear if the above equation is integrated over the hollow 
region in SST model. Eq. (|3.22p denotes the first law. 
Second law is satisfied as mentioned in Eq. (|3.10p . 

Third law is satisfied by definition of Tsst{x), if the third law of ordinary equilibrium thermodynamics holds for 
inner and outer black bodies. 

Furthermore, as for ordinary equilibrium thermodynamics, we can find the free energy density is concave with 
intensive state variables, 

From the above, we have obtained a consistent two-temperature steady state thermodynamics for radiation 
fields. 



4 Black hole evaporation with energy accretion 

Now we apply the steady state thermodynamics for radiation fields to the NE model. Contents of this section are 
based on llj. This section is not necessary to read next SecHJ Readers interested in the generalized second law 
may skip over this section. 
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4.1 From SST to NE model 

The NE model is obtained from the SST model by setting the system spherically symmetric and assigning Eq. (|2.ip 
to the inner black body as its equations of states. Then the inner black body is regarded as a black hole, and the 
net energy flow from black hole to heat bath via radiation fields causes the black hole evaporation. 

Before proceeding to the NE model, let us review here about existing works. In the framework of ordinary 
equilibrium thermodynamics, equilibrium states of black hole in a heat bath have already been investigated. It has 
already been revealed that an equilibrium state of the total system composed of a black hole and a heat bath is 
unstable for a sufficiently small black hole and stable for a sufficiently large black hole [171 HI] • If the instability 
occurs for a small black hole and the system starts to evolve towards the other stable state, there are two possibilities 
of its evolution: The first is that, due to the statistical (and/or quantum) fluctuation, the temperature of heat bath 
exceeds that of black hole and a net energy flow into black hole arises. Then the black hole swallows a part of 
heat bath and settles down to a stable equilibrium state of a larger black hole in heat bath. The second possible 
evolution is that, due to the statistical (and/or quantum) fluctuation, the temperature of heat bath becomes lower 
than that of black hole and a net energy flow from black hole arises. Then the black hole evaporates and settles 
down to some other stable state. However we do not know the detail of end state of the second possibility, since 
the final fate of black hole evaporation is an unresolved issue at present. 

When one distinguishes the phase of the equilibrium system by a criterion whether a black hole can exist stably 
in an equilibrium with a heat bath or not, the phase transition of the system occurs in varying the black hole radius. 
This phenomenon is known as the black hole phase transition [171 IH] ■ So far cosmological constant is not considered. 
But the black hole phase transition has also been found for asymptotically anti-de Sitter black holes, which is known 
as Hawking-Page phase transition [19j . However, we do not consider cosmological constant throughout this article. 

The black hole phase transition is one of interesting issues in black hole thermodynamics. However this section 
concentrates on a black hole evaporation in a heat bath after an instability of equilibrium occurs. We investigate 
a detail of black hole evaporation in the framework of NE model and try to extract some insight into the final 
fate of black hole evaporation. Readers interested also in the black hole phase transition may see [TT] in which its 
equilibrium and nonequilibrium versions are also discussed. 

4.2 Energy transport in the NE model 

We discuss energetics of NE model. Total energy of the whole system is 

Etot ■= Eg + Efi+ Erad , (4.1) 

where Eg is the energy of black hole given in Eq. (|2.ip , E^ is the energy of heat bath defined by ordinary thermo- 
dynamics, and Erad is the steady state energy of radiation fields given in Eq. (|3.9a[) . 

Erad = 4a' {Gg + Gh T^) , (4.2) 

where 

Gg -.^ [ dx'^ggix) , Gh -.^ [ dx^'ghix), (4.3) 



where gg{x) is the solid angle (divided by An) covered by directions of particles which are emitted by black hole and 
come to a point x (see figslT]and[2]), and gh{x) is defined similarly by particles emitted by heat bath. By definition 
gg{x) + gh{x) = 1 holds, and consequently Vrad Gg + Gh gives the volume of hollow region. Furthermore, since 
the black hole is concentric with the hollow region, we obtain gg{x) = ( 1 — cos 6* ) /2 and gh{x) = ( 1 + cos 6* ) /2, 
where 9 is the zenith angle which covers the black hole at a point of radial distance r (see right panel in Fig[T]). 
Then Gg and Gh are expressed as 



G, 



3 

r 

Gh 



^ r n3 n3 / d2 n2 \ 3/2 
Hi, — -ft. — I rii, — It, 



'h ^'-g V ''h ^ "-g , 



(4.4a) 
(4.4b) 



where Rh is the outermost radius of hollow region (see right panel in Fig[T|). 

To understand the energy flow in NE model, we divide the whole system into two sub-systems X and Y as 
follows: Sub-system X is composed of the black hole and the "out-going" radiation fields emitted by black hole, 
and sub-system Y is composed of the heat bath and the "in-going" radiation fields emitted by heat bath (see left 
panel in FiglT]). The sub-system X is a combined system of components of NE model which share the temperature 
Tg, and Y is that which share the temperature Th- Then the total energy is expressed as 

Etot = Ex+Ey, (4.5) 
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where Ex and Ey are respectively the energies of sub-systems X and Y, 



Ex^E. + E^fJ, , Elll.= Aa'G,T^ (4.6a) 
Ey=E, + El':l , El':l:=4a'G^T;^, (4.6b) 

where -E^od ^^"^ -^rod respectively the energies of out-going and in-going radiation fields. It is easily found that 
Ex has no T^j-dependence, while Ey has Tg- and T/j-dependence. The energy flow in NE model can be understood 
as an energy transport between sub-systems X and Y. Because this energy transport is carried by the out-going and 
in-going radiation fields, the Stefan-Boltzmann law works well to give an explicit expression of energy transport, 



dE 



-o'{Tl-n)Ag (4.7a) 
^ = a'{Tl-n)Ag, (4.7b) 

where Ag = AttR'^ is the surface area of black hole, and t is a time coordinate which corresponds to a proper time 
of a rest observer at asymptotically flat region if we can extend the NE model to a full general relativistic model. 
Because some particles emitted by heat bath are not absorbed by the black hole but return to the heat bath (see 
right panel in Fig[T]), the effective surface area through which Y exchanges energy with X is equal to the surface 
area of black hole. Therefore Ag appears in Eq. (|4.7bp . Furthermore Eas. (l4.7p are formulated to be consistent with 
the isolated setting of the NE model, Etot = constant. 

It is useful to rewrite the energy transport (j4.7p to a more convenient form for later discussions. By Eas. (|2.2p . 
and (|46l) . we find Eq. (|4jp becomes 

Cx^=-J , CxC^'^^^{Cx + C^y'^) J, (4.8) 
J:=a'{T^-T^)Ag, (4.9) 

Cx ^=C, + Ci:^, (4.10a) 



where 
and 



C'Z -^^\^a'GgTl + — [R.g-^Rl~Rl)Tg (4.10b) 



4 



C?" := '^ = C„ + 1(>o'GhTI (41l)d) 

01 h 

Ch ■■= 777;^ (>0), (4.10e) 

where Cg = — 27ri?g is given in Ea. (j2.2p and it is assumed for simplicity that Rh = constant and E^ depends on 
Th but not on Tg. Ch is the heat capacity of heat bath, and we assume Ch = constant > for simplicity. Cy ^ is 
the heat capacity of sub-system Y under the change of Tg with fixing Th, and Cy'^ is that under the change of Th 
with fixing Tg. C^^)j is the heat capacity of out-going radiation fields, and Cx is the heat capacity of sub-system X. 
In analyzing the nonlinear differential equations (|4.8p . behaviors of various heat capacities (|4.10p are used. Some 
useful properties of these heat capacities are explained in next Subsec lOl 

From the above, we find that an inequality Cx + C'y < has to hold in order to guarantee the validity of 
NE model. To understand this requirement, consider the case Tg > Th for the first. Due to the temperature 
difference, energy flows from black hole to heat bath via radiation flelds, dEg < and dEh > 0. Then dTg > and 
dTh > hold due to Cg < and Ch > 0. Recall that the whole system is isolated, Etot = constant, which means 
{Cx + Cy^) dTg + Cy^ dTh = 0. Therefore, because of Cy^ > by definition, it is concluded that the inequality 
Cx + Cy'^ < must hold. And an inequality Cx < follows immediately due to Cy ' > by definition. The 
similar discussion holds for the case Tg < Th, and gives the same inequality. Hence the following inequality must 
hold in the framework of NE model, 

Cx + C^y'^ <0 (=» Cjf <0). (4.11) 
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This inequality is the condition which guarantees the vaHdity of NE model. A more detailed property of the 
combined heat capacity Cx + Cy'^ is explained in next Subsec l4.3.4[ which shows that inequality (|4.1ip can hold 
for a sufficiently small Th. Therefore we assume Th is small enough so that the inequality (|4.1ip holds. 

Concerning the validity of NE model, what the quasi-equilibrium assumption implies is important. This as- 
sumption requires the time evolution is not so fast. Therefore, when a black hole evaporates, the shrinkage speed 
of black hole surface is less than unity, 



dt 



< 1 . 



(4.12) 



This inequality is also the condition which guarantees the validity of NE model. Hence, in the framework of NE 
model, our analysis should be restricted within the situations satisfying conditions (|4.1ip and (|4.12p . 

It is helpful for later discussions to consider what a violation of validity conditions (|4.1ip and (|4.12p denotes. 
Firstly consider if condition (|4.1ip is not satisfied. Then the system, especially the radiation fields, can never be 
described with steady state thermodynamics. The radiation fields are neither equilibrium nor steady (stationary 
nonequilibrium) . This means that the radiation fields should be in a highly nonequilibrium dynamical state. 
Furthermore the quasi-equilibrium assumption is violated, because it is this assumption that lead us to utilize the 
steady state thermodynamics. Therefore highly nonequilibrium radiation fields make a black hole dynamical, and 
the black hole can not be treated by equilibrium solutions of Einstein equation. Next consider if condition (|4.12p 
is not satisfied. Then the black hole evolves so fast that the quasi-equilibrium assumption is violated. The black 
hole can not be described by equilibrium solutions of Einstein equation, but described by some unknown dynamical 
solution. Therefore, because the source of radiation fields becomes dynamical, radiation fields evolve into a highly 
nonequilibrium dynamical state and the steady state thermodynamics is not applicable. Hence, when one of the 
conditions (|4.1ip or (|4.12p is violated, the system evolves into a highly nonequilibrium dynamical state which can 
not be treated in the framework of NE model. 

4.3 Properties of some heat capacities 

This subsection summarizes properties of various heat capacities (|4.10p which we will use in remaining subsections. 
4.3.1 Cx{Rg) as a function of Rg 

Here we show a behavior of heat capacity Cx of sub-system X as a function of black hole radius Rg. By Eqs. (|2.ip 
and (|4.4p . Cx{Rg) is rewritten into the following form. 



Cx{x)^-2i,Rix' + — f{x), 



where x := Rg/Rh and 



fix) ■■- 



4 1 
3 ^ 



l-x-" 



(i-x^f^] +i (^x~VT^) 



(4.13) 



(4.14) 



By definition, < Rg < R^, i.e., < x < 1, and we find 



Cx ^ oo 



C 



X 



-27:Rl 



87r2 



as x 



at a; = 1 



(4.15) 



By Eq. fO)) and Eq. ([^ which denotes Rh > 1, we find Cx{x = 1) < holds for the NE model. 
The differential of Cx{x) is 



dCxjx) 
dx 



-AttRI X + -— r 



a' dfix) 



where 



df{x) _ 4-3a;2-4Vl-.T2 



87r2 dx ' 
df{x 



dx 



dx 
dfjx) 

dx 



< for X < 



> for X > 



2V2 ' 



(4.16) 



(4.17) 



where 2\/2/3 ~ 0.943. A schematic graph of Cx{Rg) is shown in FigO (left panel), where Rg is the solution of 



Cx - 0, 



Rq '■— Rh X 



■.^{x\Cx{x)^Q} 



(4.18) 



Because of Cx{x = 1) < 0, equation Cx = has only one solution. 
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Cx(Rg) 



R„ ^ 0.055 (NRh) 



1/3 



Cx(Rh) = Q'|^(Rh) + constant 
3j2_ 




Ro 




Figure 3: Left panel is for Cx{Rg) as a function of Rg, where 2^2/3 ~ 0.943. Right is for Cx{Rh) as a function of 
Rh, where 3v^/4 ~ 1.06. Cx{Rg = Rh) < holds due to a' ~ 10 and I < Rg < Rh- 



Finally, we estimate the value of x. Since a; < 1 by definition, we apply the Taylor expansion, (1 — .t^)" 
1- ax^ +0{x'^),to Eq. (11331) and obtain 



Cxix) = -27ri?2 



x^-e - (l-^+0{x^] 



(4.19) 



where e = a'/Wn^Rl = iV/19207ri?2. By Eq.([Lni) and Eq.fESl) which denotes Rh > 1, we find £ < 1. Then Cx(a 



becomes 



Cx (x) = -2ttRI - (x^ -e + eOix)) 

X 



Hence, by taking leading terms of x and e, we find an approximate expression for x as 

a/3 



a; ~ e 



This gives an approximate value of Rg, 



Rn 



( NRhY^' 
V19207r 



0.055 X {NRh f^^ 



(4.20) 
(4.21) 

(4.22) 



4.3.2 Cx{Rh) as a function of Rh 

Here we show a behavior of heat capacity Cx of sub-system X as a function of outermost radius Rh of the hollow 
region. We make use of the calculations done in previous subsection. According to Eq. (|4.13p . Cx{Rh) is expressed 
as 



Cx{R; 



,,-~2nRl + cll^,iRh) , d/J^iRh)^^ 



(4.23) 



where x := Rg/Rh, < a; < 1 by definition, and f{x) is given in Ea. (|4.14p . This denotes Cx{Rh) behaves as 
c'fad^Rh) + constant. The differential becomes 



dCxjRh) 
dRh 



Stt' 



df{x) 
dx 



Using Eq. (|4.17p . we find 



dCx 
dRh 
dCx 
dRh 



a 

8^ 



CX) as ^ i?g + (a; ^ 1 + 0) 
as, Rh ^ GO (x ^ 0) 



(4.24) 



(4.25) 



Hence referring to limit values (|4.15p . a schematic graph of Cx{Rh) is obtained as shown in Figl3] (right panel). 
Cx{Rh) is monotone increasing for Rh > (3\/2/4)-Rg ~ l.OdRg. 

4.3.3 Proof of the inequality Cg/Cx > 1 

Here we prove inequality Cg/Cx > 1 under the condition Cx < (see condition (|4.1ip ). We make use of the 
calculations done in Subsec H.3.11 By Eq. (|4.13p and definition of Cx given in Ea. (|4.10p . C'fJ^ is expressed as 
Cll^ix) = {a'/8TT^) f{x), where a: := Rg/Rh, < a; < 1 by definition, and f{x) is given in Eq. (|iTl)) . Then 
Eq. gTTl) indicates C'f^^ix) > C^f^(2V2/3), and Eq.gHH gives /(2V2/3) ~ 0.845 > 0. Therefore we find C^l^ > 
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Figure 4: Numerical plots of Cx + Cy as a function of Rg with N = 100. Left three panels are for i?^ = 5 and 
Th = 0.2, 0.096 and 0.06 downwards. Right three panels are for R,, = 10^ and T,, = 0.05, 0.00075 and 0.0005 
downwards. As far as the author checked, the same behavior is observed for every value of Rh > 1. C'x + ' < 
holds for a sufficiently small T\. 



for < a; < 1, which is consistent with a naive expectation that an ordinary matter like radiation fields has a 
positive heat capacity. 

On the other hand a required condition Cx {— Cg + C^^^^) < indicates < C'^^J^ < \Cg\. This gives \Cx\ = 
\Cg\- C^^Jj^ < \Cg\. Hence we find the inequality, Cg/Cx > 1- 

4.3.4 Cx + C{?^ as a function of Rg 

Here we show a behavior of combined heat capacity Cx + Cy'' as a function of black hole radius Rg. We make use 

of the calculations done in Subsec H.3.11 By Eq. (|4.13p and definitions of Cx and Cy given in Eg. (QUI). Cx+C^'^^ 
is expressed as 

cr' r / 

Cx{Rg) + C^'\Rg) = -2TrRlx^ + — [fix)+qx'' i^x + ^T^j^ , (4.26) 

where q := (47ri?/j T/i)^, x := Rg/Rh, < a; < 1 by definition, and f{x) is given in Eq. (l4.14p . Then, using hmit 
values (|4.15p . we find 

'Cx + C^^ ^ oo as a; ^ 

( X cr' ct' ■ (4.27) 

Cx + C^'^ = -27rRl + — + —q atx = l 



15 



The first two terms in Cx + Cy at a; = 1 is negative in the framework of NE model as discussed in Eq. (|4.15p . 
Therefore, if Th is sufficiently small, then q can be small enough so that Cx + Cy^ at a; = 1 becomes negative. 
The differential of Cx+Cy'^ is complicated and not suitable for analytical utility. Instead of analytical discussion, 

we show some numerical examples of Cx + Cy^ in FigH) It is recognized with this figure that, for a sufficiently 
small T/i, the condition (|4.1ip is satisfied for a certain range of Rg. Although FigU shows only some examples, the 
same behavior is observed for every value of > 1 as far as the author checked. Therefore the validity condition 
gHI), Cx + ^ < 0, holds for a sufficiently small Th. During a semi-classical and quasi-equilibrium stage of black 
hole evaporation {Rg{t) > 1), it is reasonable to require Tg < 1 due to Eq. (l2.1|) . and also Th < Tg < 1 for black hole 
evaporation. We assume throughout this article that Th is small enough so that the equation Cx{Rg) + CY\Rg) = 
has one, two or three solutions of Rg. 

Finally we discuss what happens along a black hole evaporation for the case that the equation Cx{Rg) + 
"*(-/?£,) = has two or three solutions. For the first consider the case of three solutions, and denotes these 
solutions by ^^^o^ -^g™^ ^'^'^ ^go increasing order, R^^^ < R^^^ < R^q. If the evaporation starts with initial black 
hole radius Rg{0) in the range RgQ < Rg{0) < Rh (see for example left-center panel in FigH]), then the evaporation 

process is treated in the framework of NE model until Rg decreases to R^q . Then, when Rg reaches i?gQ , the NE 
model becomes inapplicable to the evaporation process because the validity condition ()4.11|) is violated in the range 
R'f'^ < Rg< Rgl However we can expect Rg decreases to i?^™'' even if NE model is not applicable. Then the NE 
model becomes applicable again after Rg decreases less than i?^™''. The NE model is applied to the evaporation 
process in the range RgQ < Rg < i?^™^ . 

Next consider the case that the equation Cx{,Rg) + = has two solutions. When the black hole 

evaporation starts with initial radius larger than the larger solution of Cx{Rg) + CY\Rg) = (sec for example 
right-center panel in Fig|4]), then the similar discussion given in previous paragraph holds. The NE model is not 
applicable to the evaporation process until Rg becomes smaller than the larger solution of Cx{Rg) + Cy^^g) ~ ^■ 
However after Rg becomes less than the larger solution, the NE model becomes applicable until Rg decreases to the 
smaller solution of Cx{Rg) + c'^\Rg) = 0. 

From the above we find that, for both of the cases that equation Cx{Rg) + CY\Rg) ~ has two and three 
solutions, it is the lowest solution of that equation at which the evaporation process goes out of the framework of 
NE model and proceeds to a highly nonequilibrium dynamical stage (see last paragraph in Subsec H?^ . 

4.4 Nonequilibrium effects of energy flow 
4.4.1 General aspect of the NE model 

We discuss the black hole evaporation in the framework of NE model. To analyze energy transport equations (14. 8p 
from energetic viewpoint, we consider the energy emission rate Jne by black hole. 



where Eqs. (|2.2p and (|4.8p are used. The stronger JnE: the more rapidly the mass energy of black hole Eg decreases 
along its evaporation process. The stronger emission rate Jne denotes the acceleration of black hole evaporation. 

As mentioned in Eq. (j4.10p , we assume Rh = constant for simplicity. Rh is the parameter which controls the size 
of nonequilibrium region around black hole. To understand the nonequilibrium nature of black hole evaporation, it 
is useful to compare two situations which differ only by the value of Rh with sharing the same values of the other 
parameters of NE model, i?g, Th, Ch and N. To do this comparison, we note the following three points; firstly 
Cg < by definition, secondly Cx < given in the condition (|4.1ip , and finally that | Cx \ is monotone decreasing 
function of Rh for Rh > (3-\/2/4) Rg ~ 1.06i?g under the condition Cx < (see Subsec l4".3.2p . The first and second 
points denote Cg/Cx > 0, then the third point concludes that Cg/Cx is monotone increasing function of Rh for 
Rh > (3\/2/4) Rg under the condition Cx < 0. Hence it is recognized that, for the case Rh > (3-\/2/4) Rg, the 
larger we set the nonequilibrium region, the stronger the emission rate J^e and the faster the black hole evaporation 
proceeds. Numerical examples are shown later in Subsec l4.4.4l 

The above discussion is a comparison of NE model of a certain value of Rh with that of a different value of 
Rh. In the following subsections, we compare the NE model with the other models of black hole evaporation, the 
equilibrium model used in [17' and the black hole evaporation in an empty space (a situation without heat bath 
originally considered by Hawking in Jj). 




(4.28) 
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Figure 5: Energy emission rate by black hole —dEg/dt. Left panel is Jeq for equilibrium model extrapolated to a 
black hole evaporation, Th < Tg. Center is Jne for NE model. Right is Jempty for a black hole evaporation in an 
empty space with ignoring grey body factor. White arrow in each panel expresses energy flow from black hole. 



4.4.2 Comparison with the equilibrium model used in |17] 

In the original work [17] suggesting the black hole phase transition, only the equilibrium of the system which 
consists of a black hole and a heat bath is considered. This equilibrium model is obtained from the NE model 
by setting Rh = Rg (no hollow region) and Th = Tg (equilibrium). Obviously the equilibrium model does not 
include nonequilibrium nature of black hole evaporation, since the radiation fields disappear. Exactly speaking, the 
evaporation process is not described by the "equilibrium model" . However even in the framework of equilibrium 
model, we can find the black hole evaporation occurs for a sufficiently small black holes as mentioned in SubsecHTTJ 
By extrapolating the equilibrium model to the evaporation process, we may set Th < Tg with keeping the condition 
Rg = Rh (see left panel in FigE]). Then the energy emission rate Jeq by black hole in equilibrium model is given 
by setting Rh = Rg in Jne, 

Jeq := JNE\n,=R^ = ^' {T^ - T^) Ag ^ J , (4.29) 

where J is given in Eq.diJl). We find Jne = (Cg/Cx) Jeq- Here note that Cg/Cx > 1 is shown in Subsec l4X3l 
Therefore, when the values of Rg, Th and N are the same for the NE and equilibrium models, then Jne > Jeq 
holds. This implies that the black hole evaporation in NE model proceeds faster than that in equilibrium model. 
We can recognize that the nonequilibrium effect of energy exchange between black hole and heat bath accelerates 
the black hole evaporation. 

4.4.3 Comparison with the black hole evaporation in an empty space with ignoring grey body factor 

In the original work [1] of the Hawking radiation. Hawking considered mainly a simple situation as seen in right 
panel in fig. [5l a black hole in an empty space (a situation without heat bath) with ignoring curvature scattering 
of Hawking radiation. This describes a black hole evaporation in an empty space with ignoring the so-called grey 
body factor. There is no energy accretion onto black hole in this simple situation, and time evolution is given by 
the Stefan-Boltzmann law, dEg/dt — —a'T^Ag. Usually in most of the existing works on black hole physics, the 
time scale of black hole evaporation is estimated by assuming this simple situation. 

It is interesting to compare the NE model with the black hole evaporation in an empty space with ignoring grey 
body factor. The energy emission rate Jempty by black hole in an empty space with ignoring grey body factor is 
given by the Stefan-Boltzmann law as follows, 

Jempty '■= 0"' Tg Ag , (4.30) 

where it is assumed that matter fields of Hawking radiation is the non-self-interacting massless matter fields as for 
the NE model. Then we find ^ 

Jne = ^ ^1 - Jempty ■ (4-31) 

Recall that Tg > Th holds generally for a black hole evaporation, and Cg/Cx > 1 holds in the framework of 
NE model (see Subsecg331). Then the factor (Cg/Cx) (l - T^JT^^) may be greater or less than unity. It is not 
definitely clear which of Jne or Jempty is larger than the other. 

One may naively expect that the incoming energy flow from heat bath to black hole in NE model never enhance 
the energy emission rate by black hole, and that the relation Jne > Jempty is impossible but Jne < Jempty must 
hold necessarily. It is true if the black hole heat capacity Cg is positive. However in the NE model, Cg < as 
shown in Eq. (j2.2p and a naive sense based on ordinary systems of positive heat capacity is not always true. An 
inverse sense against the naive sense may be offered; the more amount of energy is extracted from black hole by 
heat bath, the more rapidly the black hole emits its mass energy. Furthermore the energy emitted by Hawking 
radiation is absorbed by the heat bath and affects the incoming energy flow from heat bath to black hole. The 
energetic interaction (energy exchange) between black hole and heat bath determines the energy emission rate Jne- 
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When the negative heat capacity and energetic interaction are taken into account, a naively unexpected relation 
Jne > J empty IS also possible as discussed in following paragraphs: 

To analyze the energy emission rate by black hole Jne, it is useful to recall the decomposition of the whole 
system of NE model into sub-systems X and Y, as considered in Subsec l4?2l On the other hand, the black hole 
evaporation in an empty space with ignoring grey body factor is regarded as the system obtained by removing the 
sub-system Y from the NE model. This means that, from energetic viewpoint, the black hole evaporation in an 
empty space with ignoring grey body factor can be thought of as a relaxation process of the "isolated sub-system 
X" keeping Ex = constant. Therefore, for the black hole evaporation in an empty space with ignoring grey body 
factor, the energy emission rate by black hole Jempty is the energy transport just inside the sub-system X (from black 
hole to out-going radiation fields), and no energy flows out of X. However the energy transport (|4.7p in NE model 
is the energy exchange between sub-systems X and Y. This indicates that, in the NE model, the energy Ex of X is 
extracted by Y due to the temperature difference Tg > Th, and energy flows from X to Y. The energetic interaction 
(energy exchange) between X and Y makes the black hole evaporation in NE model quit different from the black 
hole evaporation in an empty space with ignoring grey body factor. This difference is recognized significantly by 
considering a limit of the energy transport (|4.7p as follows: One may expect that the energy emission by black hole 
in an empty space with ignoring grey body factor, dEg/dt — —Jempty, should be obtained from Eq. (|4.7p by the 
limit operations, T/i — > 0, E^ — > (remove the sub-system Y) and — > oo (infinitely large volume of out-going 
radiation fields). However these operations transform Eq. (|4.7p into the set of equations, dEg/dt = — Jempty and 
= Jempty This givcs an unphysical result Eg = constant (— oo) which contradicts the "evaporation" , dEg/dt < 0. 
The black hole evaporation in an empty space with ignoring grey body factor can not be described as some limit 
situation of the NE model. 

In addition to the naive expectation Jne < Jempty, the opposite relation Jne > Jempty may be expected due 
to the negative heat capacity of black hole ()2.2|) and energetic interaction as follows: In the NE model, because the 
energy extraction {dEx < 0) occurs along the black hole evaporation {dTg > 0) , the heat capacity of X is always 

negative Cx — dEx /dTg — Cg + C^^); < 0, where C^^)j > as indicated in Subsec l4".3.3l Furthermore the larger 
the volume of hollow region, the larger the heat capacity C^^J^ and the smaller the absolute value \Cx\ because of 
Cg < 0. This implies that the more thick the hollow region, the more accelerated the increase of Tg due to the 
relation dTg = \dEx /Cx\- Therefore, for sufficiently large Rh, the energy extraction from X by Y (the increase 
of Tg) can dominate over the in-coming energy flow onto black hole. This means the energy emission rate Jne is 
enhanced by the energy extraction from X by Y, then Jne > Jempty is implied. 

The above discussion can also be supported by the following rough analysis: When a black hole evaporates in 
NE model, the temperature difference 6T :— Tg — Th should grow infinitely, 5T — > oo, due to the negative heat 
capacity of black hole. Then, because of Ea. (|4.3ip together with the facts 1 — (Tji/Tg)^ 1 (as ST oo) and 
Cg/Cx > 1 (see Subsec l4".3.3l) . the larger the temperature difference ST, the larger the ratio Jne/ Jempty Hence 
for the black hole evaporation in NE model, it is expected that the relation Jne > Jempty comes to be satisfied 
during the evaporation process even if the relation Jne < Jempty holds at initial time. Furthermore, if the relation 
Jne > Jempty holds for a sufficiently long time during the evaporation process, the evaporation time scale in NE 
model can be shorter than that in an empty space. Hence it is possible that the black hole evaporation in NE 
model proceeds faster than that in an empty space, where black holes of the same initial mass are considered in 
both cases. In next subsection, numerical examples support this discussion. 

4.4.4 Numerical example 

We show numerical solutions Tg{t) and Tfi{t) of energy transport equations (|4.8p . The initial conditions are 

Rg{0) = 100 , Th{0) = 0.0001 . (4.32) 

Rg{0) = 100 gives Tg{0) ~ 0.00079. The other parameters are set 

Ch = 1000 , TV = 100 , (4.33) 

where see Eq. (|1.3|) for N. Furthermore we have to specify the outermost radius Rh of hollow region. As mentioned 
in Eq. ()4.28|) . by comparison of a numerical solution of energy transport (I4.8|l of a certain value of Rh with that of 
a different value of Rh, we can observe the nonequilibrium effect of energy exchange between black hole and heat 
bath. The numerical results are shown in Fig[6l and the value of Rh is attached in each panel. Time coordinate r 
in this figure is a time normalized as 

r := — ^ , (4.34) 

'^empty 
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Figure 6: Numerical solutions of energy transport (|4.8p for Ch = 10^, N = 100 and the initial conditions, Rg{0) = 
100 and Th{0) = 0.0001. Horizontal line denotes the normalized time r := t/tempty- The outermost radius of hollow 
region Rh determines the size of nonequilibrium region. 



where tempty is the evaporation time (life time) of a black hole in an empty space with ignoring grey body factor. 
tempty is determined by the energy emission rate by black hole in an empty space, 

'-'empty : y^.oo j 

where no energy accretion exists due to the absence of heat bath and ignoring grey body factor, Eg corresponds to 
the mass energy of black hole and Jempty is given in Eq. (|4.30|) . This and Eq. (|2.ip give 

/ Nt \^^^ 

and we obtain 

tempty := Rg{0)^ ^ 4.02124 X 10^ , (4.37) 

where conditions (j4.32p and (|4.33p are used in the second equality. This tempty is usually adopted as the time scale 
of black hole evaporation in many existing works on black hole physics. 
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Figure 7: JNE/Jempty for each case in Figl6l The graph in left-center panel is continuous and smooth. 



Furthermore we consider the other time t]\[E a.t which one of the vaUdity conditions of NE model (|4.1ip or (|4.12p 
breaks down, 

ti,t2 Cx{h) + ci'\h) = o,v{t2) 



tNE '■= min 



1 



(4.38) 



where Cx + Cy'^ is regarded as a function of time t through Rg{t), and v{t) :— \dRg{t)/dt\ is the shrinkage speed 
of black hole radius. Each panel in FiglH] shows time evolutions of Tg{t) and Th{t) ioi < t < t^E- The black hole 
temperature at ijvB is denoted by T*, 

TgitNE)=:T:^-l-, (4.39) 
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which is also attached in each panel in FiglHl The time t^E and the other quantities obtained from our numerical 
results are listed; 



Rh/Rg{0) 


40.7426 


50 


240.758 


10* 


10^ 


E /iempty 


557.370 


1.64705 


1.00000 


0.915560 


0.536364 


Rq 


4.07277 


4.36041 


7.36293 


25.4986 


54.9354 


Jne/ Jempty{tNE) 


2001 


2294 


6540 


78435 


44029 


Rh/tNE 


1.8 X lO"'^ 


7.5 X 10-^ 


6.0 X lO-'' 


2.7 X 10-^ 


4.6 X 10-1 



(4.40) 



For cases of Rh/Rg{0) = 40.7462, 50, 240.758 and 10^, the numerical plots stopped by v{tNE) — 1, but for a case of 

Rh/Rg{0) = 10^, it stopped by Cx{tNE) + H^W-e) = 0- The second line in this list for tNE/tempty supports the 
discussion given in Subsec H.4.11 that the larger the radius Rh, the faster the black hole evaporation proceeds and 
the shorter the time Ine- The third and fourth lines in list (|4.40p for R* and Jne/ Jempty{iNE) give an important 
information in next subsection. The lowest line in list (j4.40p for Rh/tNE shows that our numerical results are 
consistent with the fast propagation assumption. To understand this, consider a typical time scale trad in which one 
particle of radiation fields travels in the hollow region from black hole to heat bath. This trad is given by trad ■= Rh, 
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since the radiation fields are massless. If trad < tNE, then the fast propagation assumption is reasonable. In fact, 
the ratio trad/ tNE (= Rh/tNE) shown at the lowest line indicates the validity of fast propagation assumption. 

Concerning a case Rh = 20 Rg{0), it is helpful to recognize that our conditions (|4.32p and (|4.33p give 2 tt x 100^ > 
1000 which denotes \Cg\ > Ch- This \Cg\ > Ch is the condition for stable equilibrium of black hole and heat bath 
in the framework of the equilibrium model (see [17j or [TT]). Therefore, if the nonequilibrium region is ignored and 
the equilibrium model used in [T7] is considered with the same setting parameters of Eas. (|4.32p and (|4.33p . then 
the black hole stabilizes with heat bath to settle down into a total equilibrium state Tg — 0. Indeed, for 

the case Rh = 20i?g(0), the hollow region is not large enough and the black hole stabilizes with heat bath. Hence 
the occurrence of accelerated increase of temperature Tg in FiglH] is obviously the nonequilibrium effect of energy 
exchange between black hole and heat bath. 

Then as the outermost radius Rh of hollow region is set larger and the nonequilibrium region becomes larger, 
the black hole evaporation comes to be observed as an accelerated increase of temperature difference Tg — Th for 
the case Rh = 40.7426 i?g(0). Ine is very longer than tempty in this case. Furthermore as Rh is set larger, the 
time tNE becomes shorter and we find tNE — tempty for Rh — 240.758 i?g(0). At last tNE becomes shorter to 
tNE — 0.536364: tempty for Rh — 10^i?g(0). If we set Rh > 10^, the combined heat capacity becomes positive 
Cx + C'y ■* > which violates the validity condition (|4.1ip . A black hole with a large nonequilibrium region of 
Rh > 10^ can not be treated in the framework of NE model, and, as discussed in the last paragraph in Subscc l42l 
a black hole evaporation for such case is a highly nonequilibrium dynamical process in which the black hole can not 
be treated with equilibrium black hole solutions of Einstein equation. The larger the nonequilibrium region, the 
faster the black hole evaporation process evolves into a highly nonequilibrium dynamical stage. 

Fig[7] shows time evolutions of the ratio of energy emission rates Jne/ Jempty given by Eq. (|4.3ip . This figure 
indicates that, when a black hole evaporates, the relation Jne > Jempty comes to hold even if the converse relation 
Jne < Jempty holds at initial time. Therefore the discussion given in previous Subsec l4.4.3l is supported. For 
the case Rh = 20Rg{0) in which a black hole stabilizes with heat bath, the energy emission rate Jne disappears, 
Jne/ Jempty 0, as easily expected by the behavior Tg Th Jne 0. 



4.5 Beyond the NE model I: abrupt catastrophic evaporation 

When the nonequilibrium region around black hole is not so large, the evaporating black hole is well approximated 
to equilibrium solutions of Einstein equation and the NE model is applicable to such quasi-equilibrium evaporation 
stage. However the condition (|4.11[) or (|4.12p is violated at the time tNE and the evaporation process becomes 
highly nonequilibrium dynamical stage (see last paragraph in Subsec l4?2)) . Therefore, since a relation i?* > 1 is 
expected from the list (|4.40p . we find the semi-classical evaporation stage (Rg > 1) is divided into two stages, 
quasi-equilibrium one and highly nonequilibrium dynamical one. The former is described by the NE model, but the 
latter is beyond the range of NE model. This subsection extrapolates NE model to the latter stage and suggests a 
specific nonequilibrium phenomenon. 

On the other hand, for the equilibrium model used in 17J and the black hole evaporation in an empty space 
with ignoring grey body factor, the highly nonequilibrium dynamical stage does not exist and the semi-classical 
stage is always described as the quasi-equilibrium stage. This is explained in next Subsec H.5.11 Then the highly 
nonequilibrium dynamical stage and a suggestion about that stage in the framework of NE model are discussed in 
SubsecgXll 



4.5.1 On semi-classical evaporation stage, except for NE model 

This subsection shows that the highly nonequilibrium dynamical stage does not occur for the equilibrium model 
used in |17j and for the black hole evaporation in an empty space with ignoring grey body factor. 

For the first consider the equilibrium model used in [17]. As explained in Subsec l4.4.21 this model is obtained 
from NE model by setting Rh = Rg and Tg = Th- Exactly speaking, the evaporation process is not described by 
the "equilibrium model" . However by extrapolating the equilibrium model to the evaporation process, we may set 
Th < Tg with keeping the condition Rg — Rh- Then the energy transport equations from black hole to heat bath 
are given by 

^ = -a'{T^-T;t)Ag , ^=a'{T^-T;^)Ag, (4.41) 
where Ag = AnRg is the surface area of black hole. The equation of dEg/dt together with Eq. (l2.ip give 

= 2a'{T^-T,t)Ag. (4.42) 

Since equations of states (|2.ip are used, the quasi-equilibrium assumption is also necessary here and ti < 1 is required. 
The inequality v < 1 corresponds to the validity condition (|4.12p of NE model. Obviously v = 1 occurs for non-zero 



dRg 


= 2 


dEg 


dt 




dt 
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radius Rg > 0. On the other hand, from Eq. (l2.2|) . vanishing heat capacity of black hole Cg = corresponds to 
zero radius Rg = 0. This means that the equilibrium model has no validity condition which corresponds to the 
condition (|4.1ip of NE model. Hence we can recognize that, if w = 1 occurs for a black hole of semi-classical size 
Rg > 1, it is concluded that the highly nonequilibrium dynamical stage occurs at semi-classical level. But if ?; = 1 
does not occur for a semi-classical black hole, then the highly nonequilibrium dynamical stage does not occur at 
semi-classical level in the framework of equilibrium model. 
The validity condition u < 1 is rewritten as 

a'(T4-r,^) <2^r2, (4.43) 

and this gives 

^'>8^(^ + Vl + (2Wf) (4-44) 

where f3 := 27r/cr' = 240/7riV and Rg = l/'iirTg is used. Here Eq.lLH]) gives /3 ~ 1. Recall that Tg > Th holds 
generally for any evaporation process and Tg < 1 holds due to the quasi-equilibrium assumption. Then we find 
Th <l and (2r^//3)^ < 1. Therefore we can approximate inequality (|4.44[) to Rg > l/(167r2/?) ~ 10"^, where /3 ~ 1 
is used. This gives 

Rg > 0.1 . (4.45) 

This denotes that the fast evaporation v = 1 occurs at Rg ~ 0.1. Hence, as discussed in previous paragraph, the 
highly nonequilibrium dynamical stage does not occur in the framework of equilibrium model used in |17j . 

Next consider a black hole evaporation in an empty space with ignoring grey body factor. The Stefan-Boltzmann 
law gives Ea. (|4.35p . This together with equations of states (j2.ip give 

= 2<y'T^Ag . (4.46) 

Hence, following a similar discussion given in previous paragraph with setting Th — 0, we require v < 1 and obtain 
Rl > l/(167r2/3) ~ 10-2, where /3 ~ 1 is used. This gives 

Rg > 0.1 . (4.47) 

This denotes that the fast evaporation v = I occurs at Rg ~ 0.1, and that the highly nonequilibrium dynamical 
stage does not occur for a black hole evaporation in an empty space with ignoring grey body factor. 
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4.5.2 Abrupt catastrophic evaporation 

Let us point out again a numerical evidence shown in list (|4.40p that the black hole radius R* at time t^E is 
greater than unity R* > 1. According to the discussion in last paragraph in SubsecH^l we can expect a highly 
nonequilibrium dynamical stage of evaporation process will occur at semi-classical level R* > 1 in the framework 
of NE model. After the time t^E, the black hole and radiation fields should be described as highly nonequilibrium 
dynamical ones. 

In the following discussion, we make two steps: Firstly, to confirm the numerical evidence, we show R* > 1 
analytically. Secondly, a physical implication of _R* > 1 is discussed and a specific nonequilibrium phenomenon is 
suggested. 

For the first, we analyze the energy transport equations (|4.8p . Due to definition (|4.38p of time t^E, we consider 
two cases, Ine = ti and Ine = ^2, where ti and ^2 are given in definition (|4.38p . But before proceeding to 
the analysis of these cases, we should point out the following: As explained at the end of Subsec l4.3.4[ equation 
Cx{Rg) + d^XRa) = has one, two or three solutions of Rg for a sufhciently small Th- However even if there 
are two or three solutions of Rg for our choice of Th, it is the lowest solution at which a highly nonequilibrium 
dynamical evaporation stage starts towards a quantum evaporation stage. Hence the time ti in definition (j4.38p is 
the lowest solution of equation Cx{t) + Cy''(0 = 0. 

Here we estimate the order of i?*. Consider the case Ine = ti, where R* = Rgih) and Cx{R*g) +C'y ''(-Rg) = 0. 

Because of Cy^ > by definition, Cx{R*g) = — C'y^(^g) < holds. Consequently, according to a behavior of 

Cx{Rg) explained in SubsecgXT] (left panel in FigED, we find R*g> Rg:^ 0.055 x {NRhf^^. Hence together with 

Eq. (|1.3p . we find -Rg > 1 for the situation Rh ^ 60. And next consider the other case t^E — ^2, where \Rg{t2)\ ~ 1. 

Because of t2 < ii, we find CxiH) + C'y ■'(t2) < 0. Then, because it is assumed that Th is small enough so that 
the vahdity condition (14. Up holds, we find by Subsec HXil (FigH)) that i?* = Rg{t2) > Rgo holds, where Rgo is the 

lowest solution of Cx{Rga) + CY\Rgo) = 0. Therefore, following the same discussion given for the case t^E = ^i, 
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we obtain R* > Rgo > i?g ~ 1 for the situation Rh > 60. In summary, the black hole radius R* at time Ine 
is greater than unity R* > 1, when the black hole evaporates in the framework of NE model under the condition 
Rh ^ 60. Here we have to note two points: First is that the condition Rh > 60 is not a necessary condition but a 
sufficient condition for i?g ^ 1, and there may remain a possibility that Rg^l holds even if Rh < 60. Second point 
is that, if the nonequilibrium nature of black hole evaporation is not taken into account, the radius R* can not be 
greater than unity but it becomes less than Planck length as seen in Eas. (|4.45p and (|4.47p . The relation i?* > 1 is 
a peculiar property of the NE model. 

We proceed to the second part of this subsection, an implication of the above result, R* > 1. Recall a highly 
nonequilibrium dynamical stage of evaporation process occurs after the time t^E- Because of i?* > 1, a semi- 
classical (but not quasi-equilibrium) discussion is available for the highly nonequilibrium dynamical stage while 
the black hole radius shrinks from R* to Planck length Ipi := 1. Then it is appropriate to consider that the mass 
energy of black hole evolves from E* {— Rg/2) to Ep := lpi/2. Energy difference AEg := E* — Ep is emitted 
during highly nonequilibrium dynamical stage. Furthermore, for example, Fig[7]and fourth line in list (|4.40p of our 
numerical example imply a very strong luminosity Jn e of Hawking radiation in the NE model in comparison with 
the luminosity Jempty in the evaporation in an empty space with ignoring grey body factor. On the other hand 
Jempty is vcry strong as explained in §1 of [T]. Hence, Jne may be a huge luminosity. The energy emission by a 
black hole in NE model may be understand as a strong "burst" . 

In addition to the luminosity of Hawking radiation, we consider the duration dtdyn of the highly nonequilibrium 
dynamical stage. Since the shrinkage speed of black hole radius v := \dRg/dt\ is approximately unity w ~ 1 during 
highly nonequilibrium dynamical stage (see condition (|4.12p ). the duration is estimated as St^yn ^ Rg/v ^ R*, and 
the following relation is obtained, 

^ i < M = ^ ' « 0.025, (4.48) 

tempty ^empty tempty 12807r_/?(^ (0) AORg{Oj 

where Ea. p.3p is used, and, since the initial radius Rg{0) should be large enough to consider a semi-classical 
evaporation stage, we introduced relations R* < Rg{0) and 1 ^ Rg{0)- This denotes dtdyn ^ tempty Furthermore, 
for example, we find the shortest Ime = 0(0. 1) x tempty from list (|4.40p . This together with (I4.48P imply 5tdyn ^ 
ijv_B- Hence it seems reasonable to consider that Stdyn is very shorter than tjvB- (For example it seems that 
the tangent dTg/dt seen in each panel in FigE] becomes very large quickly as t — > ijvB: and Tg will reach Planck 
temperature quickly just after t^vB-) Hence it is suggested that the energy AEg := E* — Ep bursts out of black 
hole with a very strong luminosity within 6tdyn which is negligibly short in comparison with t^E- 

From the above, we suggest the following: When a black hole evaporates in the framework of NE model under 
the condition Rh ^ 60, a quasi-equilibrium evaporation stage continues until tNE- Then a highly nonequilibrium 
dynamical evaporation stage occurs at tNE- In that stage, a semi-classical black hole of radius -Kg > 1 evaporates 
abruptly (within a negligibly short time scale dtdyn) to become a quantum one. This abrupt evaporation in the 
highly nonequilibrium dynamical stage is accompanied by a burst of energy AEg. We call this phenomenon the 
abrupt catastrophic evaporation at semi-classical level -Rg ^ 1, where "catastrophic" means the shrinkage speed 
of black hole radius is very high v ^ 1 and the energy AEg bursts out of black hole with a huge luminosity 
Jne ^ Jempty within a negligibly short time scale 6tdyn tNE- 

The above discussion is based on the NE model. As shown in previous subsection, for the equilibrium model used 
in [1] and the black hole evaporation in an empty space with ignoring grey body factor, the black hole radius becomes 
Planck size before the shrinkage speed of black hole radius reaches unity. The highly nonequilibrium dynamical 
stage and the abrupt catastrophic evaporation at semi-classical level do not occur in those models. Hence the abrupt 
catastrophic evaporation at semi-classical level seems to be a specific nonequilibrium phenomenon suggested by NE 
model. 

Here we discuss about a black hole evaporation in an empty space in a full general relativistic framework. 
Note that, even if a black hole is in an empty space, there should exist an incoming energy flow onto the black 
hole due to the curvature scattering. When the curvature scattering is taken into account for the case of a black 
hole evaporation in an empty space, we can interpret the whole system as if a black hole is surrounded by some 
nonequilibrium matter fields which possess outgoing and incoming energy flows of Hawking radiation under the 
effects of curvature scattering. Furthermore, since the curvature scattering occurs whole over the spacetime, it is 
expected that the nonequilibrium region is so large that a condition corresponding to Rh ^ 60 in NE model holds. 
Hence, if the NE model is extended to a full general relativistic model, we can expect that a black hole evaporation 
in an empty space can be treated in the framework of full general relativistic version of NE model (with removing 
the heat bath), and that the abrupt catastrophic evaporation at semi-classical level may occur as well since the 
nonequilibrium region is sufficiently large. 

Finally we estimate a typical time scale of black hole evaporation with energy accretion. It is reasonable to 
consider the duration of quantum evaporation stage is about one Planck time. Then, the time scale of black hole 
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evaporation tev is estimated as 



t, 



ev — 



tNE + Stdyn + 1 ~ tj\[E ■ 



(4.49) 



The time t^E gives a typical time scale of black hole evaporation with energy accretion. 

4.6 Beyond the NE model II: final fate of quantum black hole evaporation 

So far we have considered semi-classical evaporation stages in the framework of NE model, and found it consists of 
two stages, quasi-equilibrium one and highly nonequilibrium dynamical one. This subsection discusses the quantum 
evaporation stage following the highly nonequilibrium dynamical stage. 

Concerning quantum black hole evaporation, the so-called information loss paradox is an interesting and im- 
portant issue [inj: If the black hole mass is radiated out completely by the Hawking radiation, then the thermal 
spectrum of Hawking radiation implies that only a matter field which is in a thermal equilibrium state may be 
left after the evaporation. This implies that the initial condition of black hole formation (gravitational collapse) is 
completely smeared out. For example, even if the initial state of collapsing matter is a pure quantum state, the 
final state after a black hole evaporation must be transformed to a thermal state. This contradicts the unitary 
invariance of quantum theory. This is the information loss paradox. 

The study on the final fate of black hole evaporation is usually carried out in the context of quantum gravity 
theory. However, since no complete theory of quantum gravity has yet been constructed, it is meaningful to some 
extent to study the final fate of black hole evaporation with an appropriate model of black hole evaporation without 
referring to present incomplete quantum gravity theories. Therefore, we use the NE model and try to extract what 
we can suggest about the final fate of black hole evaporation. 

In this subsection, we consider whether some remnant remains after the quantum evaporation or not. If a remnant 
remains, then it is implied that the information loss paradox does not exit, because the complete evaporation may not 
be true and the remnant may preserve the information about initial condition of black hole formation to guarantee 
the unitary evolution of the system. For the time being, we assume that a black hole evaporates out completely and 
only equilibrium radiation fields remain at the end state of quantum evaporation stage. If a contradiction results 
from this assumption, we may conclude that a remnant will remain after a black hole evaporation. This subsection 
aims to suggest a necessity of some remnant by the reductive absurdity. 

As discussed at the end of previous subsection, it is expected that the abrupt catastrophic evaporation at semi- 
classical level occurs in a full general relativistic framework. Therefore we consider the case Rh > 60 in the NE 
model which denotes the occurrence of abrupt catastrophic evaporation. Then, under the assumption of complete 
evaporation of black hole, we can draw a scenario in the framework of NE model as follows: 

As a quasi-equilibrium stage of black hole evaporation proceeds until t^E, the temperature Tg approaches a 
critical value T* = l/iirR*, where R* := Rg{tME)- Then a highly nonequilibrium dynamical stage and a quantum 
evaporation stage follow successively. Under the assumption of complete evaporation, these successive stages are 
described as one phenomenon; a black hole emits completely its mass energy E* {— R*/2) within a negligibly 
short time scale Stdyn + 1 ^tvb), and equilibrium radiation fields remain. Hereafter in this subsection until a 
contradiction will be derived, the abrupt catastrophic evaporation under the assumption of complete evaporation 
means those successive stages, highly nonequilibrium dynamical one and quantum evaporation one. Then the 
abrupt catastrophic evaporation is described by a replacement of a black hole of radius R* by equilibrium radiation 
fields of volume V* = (47r/3) i?*'^. In this replacement, thermodynamic states of heat bath and nonequilibrium 
radiation fields around the region of V* are not changed (see the upper part in Fig[8]). Here we assume the energy 
conservation so that the energy of equilibrium radiation fields in volume V* equals the mass energy E* of black 



When a sufficiently long time has passed after the abrupt catastrophic evaporation (under the assumption of 
complete evaporation), the whole system reaches an equilibrium state in which a heat bath and radiation fields in 
the hollow region have the same equilibrium temperature. However, without considering such totally equilibrium 
end state of the whole system, but with considering the states just before and just after the abrupt catastrophic 
evaporation, we can discuss whether a remnant remains or not after a quantum black hole evaporation. 

It is reasonable to require the increase of total entropy along black hole evaporation, because of the isolated 
condition of the whole system in NE model. Therefore S^^^ < S^Jf holds, where S^^f is the total entropy of the 
whole system at time t^E just before the abrupt catastrophic evaporation, and S^J^ is the total entropy of the whole 
system just after the complete evaporation of black hole. Because the abrupt catastrophic evaporation is described 



hole, 




(4.50) 
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Figure 8: If a black hole evaporates out completely, when a black hole temperature Tg reaches T*, the black hole is 
suddenly replaced by radiation fields of temperature T"*^^ which is determined by the energy conservation. However 
it is impossible from entropic viewpoint, and a remnant should remain. But its equations of states are unknown. 




Figure 9: Schematic graph of AS'j*^^. A prohibition against complete evaporation of black hole is suggested. 



by a simple replacement of a black hole by equilibrium radiation fields without changing thermodynamics states of 
heat bath and nonequilibrium radiation fields, the entropy difference is given by 

^'S'tot •= ^tot ~ ^tot — ^rad ~ i (4-52) 

where S* is black hole entropy of temperature T* and S*^^ is equilibrium entropy of radiation fields of volume V* 
and temperature T*^^. If the complete evaporation of black hole is true of the case, a relation AS'j*^^ > must hold. 
However it is impossible as follows: 

Using equations of states (12. 1|) and the equilibrium entropy of radiation fields (16cr'/3) T*j^ V^* together with 
Eq. (|4.5ip . the entropy difference becomes 

ASl,^li^^y\l'^'-nR;\ (4.53) 

A schematic graph of AS;^^ is shown in FigH where i?*o (8/27%/5^)%/]V ~ 0.075V7v is given by AS^^^ = and 
Eq.dTll). We find AS^^t > for i?; < i?;^ and AS;^t < for R*g > i?;^,. Here recall that, because of C^^^ > by 
definition and a behavior of Cx{Rg) explained in Subsec H.3.11 (left panel in Fig[3]), an inequality R* > Rg holds 
due to the validity condition (|4.1ip . where Rg is given by Cx{Rg) = 0. Furthermore, using Eq. (l4.22p . we find 
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RglR*gQ ~ (27\/5^/32(307r)i/3) x [RI/NY'^ ~ 0.734 x {Rl/Nf'^. Then, by Eq.dO]) and requirement Rh > 60 of 
the occurrence of abrupt catastrophic evaporation, we obtain 

1^>^-1- (4.54) 

This indicates AS'j*^^ < 0. Consequently, by the reductive absurdity as explained in fourth paragraph in this 
subsection, a complete evaporation of black hole is impossible (see Fig 15]). 

From the above, we give a suggestion as follows: a complete evaporation of black hole is prohibited and a remnant 
should remain after a black hole evaporation. The entropy of remnant should guarantee the increase of total entropy. 
This implies disappearance of the information loss paradox. 

Here we comment about the mass of remnant. One may think that the mass energy (or internal energy) 
of remnant may be extracted by some energetics. One may obtain a minimum entropy »S'^™^ of remnant as 
^min'^ = Sg, and a mass energy of remnant may be obtained from S^l'^\ However, because equations of states of 
remnant is unknown, it is impossible to obtain an energy of remnant from S^l^\ On the other hand, as discussed 
in Subsec HTSl (in Eq. (|4.48p ). the mass of quantum black hole is expected to be about one Planck energy. Therefore, 
as long as the NE model is extrapolated over its validity, the mass of remnant seems to be of the order of one Planck 
energy. 

In the rest of this subsection, we discuss more about the validity of discussion given in deriving inequality 
()4.54p . There are three points we are going to discuss here. For the first, we discuss about the usage of NE 
model in this subsection. If we want to analyze details of the final stage of black hole evaporation process, then 
a quantum gravity is necessary, since a black hole becomes quantum, Rg < 1. On the other hand, because of the 
quasi-equilibrium assumption, the black hole radius has to be restricted as (_Rg(0) >)Rg ^ 1 in the framework of 
NE model. However this restriction R* ^ I does not mean an impossibility of NE model for considering the final 
fate of black hole evaporation. Although the NE model is based on a classical Schwarzschild black hole, once one 
assumes a complete evaporation, then it is true that radiation fields remain and their entropy can be counted after 
a quantum evaporation process ended. This means, even though we never refer to any detail of present incomplete 
quantum gravity theories, we can compare an entropy of black hole before the start of quantum evaporation process 
with an entropy of radiation fields after the end of quantum evaporation process. Hence our analysis done above 
seems to be reasonable to approach the final fate of black hole evaporation, and does not include the "uncertainty" 
of present incomplete quantum gravity theories. Our result seems more universal than the other discussions based 
on present incomplete quantum gravity theories. 

For the second point, we turn our discussion to the starting assumption that a black hole evaporates out 
completely and equilibrium radiation fields remain. As seen above, Ea. (|4.53p has led a contradiction to deny 
the assumption and to result in necessity of a remnant. Here recall that Eq. (|4.53|) is calculated with using the 
equilibrium entropy of radiation fields. Then, one may think it is more general to modify the starting assumption 
so that radiation fields after complete evaporation are not necessarily in an equilibrium. However, as explained in 
detail in |,13 , the nonequilibrium entropy of radiation fields should be smaller than the equilibrium one 5**^^^^, since 
the equilibrium state is the maximum entropy state. Hence it is reasonable to expect that an inequality AS^^^ < 
holds stronger for a modified starting assumption which requires nonequilibrium radiation fields after complete 
evaporation of black hole. 

For the third point, we comment about black hole evaporation in an empty space (a situation without heat bath). 
As discussed at the end of Subsec HTSl the abrupt catastrophic evaporation at semi-classical level is expected to 
occur for a black hole evaporation in an empty space in a full general relativistic framework. Recall that calculations 
of entropy difference AS'j*^^ in Eq. (|4.53p do not depend on the outside of black hole, but depend only on the black 
hole entropy S* and the equilibrium entropy S*^^ of radiation fields. This means, without respect to the degree 
of nonequilibrium nature of the environment around black hole, once the abrupt catastrophic evaporation of black 
hole occurs, the same discussion as given in this subsection seems applicable to a black hole evaporation in an empty 
space. Therefore, a full general relativistic treatment of black hole evaporation in an empty space may result in 
necessity of a remnant at the end state of black hole evaporation. 



4.7 Summary of this section 

In this section, using the NE model, we found foUowings: 

• The semi-classical evaporation stage consists of two stages, quasi-equilibrium one {t < t^E) and highly 
nonequilibrium dynamical one (tjvB < t). The former is treated by the NE model, while the latter is not. If 
the steady state thermodynamics for radiation fields were not applied, the latter stage did not found. The 
existence of the latter stage is the very result of NE model. 
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• In the quasi-equilibrium evaporation stage, if the thickness of the hollow region is thin enough, then Jn e / Jempty < 
1 holds during the black hole evaporation process and t]\[E > tempty is obtained. However if the hollow region 

is thick enough, then Jne / Jempty > 1 holds during the evaporation process and t^E < tempty is obtained. The 
nonequilibrium effect of energy exchange between black hole and its outside environment tends to accelerate 
the evaporation process. Because of the negative heat capacity of black hole, an inverse sense against a naive 
sense based on positive heat capacity is offered; the more amount of energy is extracted from black hole by 
heat bath, the more rapidly the black hole emits its mass energy. 

• Duration of the highly nonequilibrium dynamical stage 5tdyn is negligibly shorter than e ■ Energy emission 
rate by black hole in NE model Jne at t^E is very stronger than that in an empty space Jempty These imply 
a huge energy burst during highly nonequilibrium dynamical stage. Such huge burst is the very result of NE 
model, and we call it the abrupt catastrophic evaporation. 

• Since the duration of quantum evaporation stage following the highly nonequilibrium dynamical stage will be 
about one Planck time, the time scale of black hole evaporation is estimated as igu — tNE + ^tdyn + 1- 
This gives t^v t^E- 

• By extrapolating the NE model to the highly nonequilibrium dynamical evaporation stage and quantum 
evaporation stage, a complete evaporation of black hole after the quantum evaporation stage is prohibited. 
This denotes a remnant of Planck size may remain at the end of quantum evaporation stage in order to 
guarantee the increase of total entropy along the whole process of evaporation. This implies a disappearance 
of the information loss paradox due to the nonequilibrium effect of energy exchange between black hole and 
its outside environment. This suggestion does not depend on present incomplete theories of quantum gravity. 

5 Physical essence of the generahzed second law 

Now we apply the steady state thermodynamics for radiation fields to the NE model, and modify it to reveal the 
physical essence which guarantees the validity of generalized second law of black hole thermodynamics. Contents of 
this section are based on [51], and written without referring to previous SecHl Readers interested in the generalized 
second law may skip over previous section. 

5.1 GSL in the context of black hole evaporation 

According to the black hole thermodynamics [l][2l[3], a classical size black hole {Rg > 1) is regarded as an object 
in thermal equilibrium, whose equations of states are given in Eq. ()2.1|) . In a general relativistic framework, the 
statement of generalized second law in the context of black hole evaporation is as follows [3] : 

Generalized second law (GSL): When a black hole evaporates, its mass energy Eg decreases and consequently 
the black hole entropy Sg decreases because of the equations of states (|2.ip . However the total entropy of 
the whole system which consists of the evaporating black hole, the Hawking radiation and any other matters 
around black hole, must increase. 

In this statement we consider all matter fields on black hole spacetime plus the black hole itself. This denotes the 
whole system under consideration is isolated. Hence the GSL requires that the time evolution of that system is a 
relaxation process and the total entropy increases. 

The original form of GSL was a conjecture that the horizon area (divided by 4) is regarded as a true black hole 
entropy |3]. However refer to the work [10] which were written after the suggestion of original GSL and revealed 
the horizon area (divided by 4) is the "equilibrium" entropy of whole gravitational field on black hole spacetime. 
Therefore, throughout this article, we consider the horizon area is a true equilibrium entropy of black hole and 
the GSL has already been proven [TUl [25] . The equilibrium black hole entropy given by the horizon area is Sg in 
Ea. (|2.ip . which is called the Bekenstein-Hawking entropy. This section aims not to prove GSL, but to reveal a 
physical essence which guarantees the validity of GSL. 

When we consider GSL in the context of black hole evaporation, the nonequilibrium nature in the outside 
environment around black hole should be taken into account. This means that, while a classical size evaporating 
black hole is described by an equilibrium solution of Einstein equation under the quasi-equilibrium assumption, the 
outside environment should be treated as a nonequilibrium matter. Then we have to define the total entropy of 
the whole system, which consists of the evaporating black hole and the nonequilibrium matter fields around black 
hole including Hawking radiation. It has already been shown in [53] that, for the equilibrium state of a black hole 
with general "self-interacting" matter fields (a heat bath), the total "equilibrium" entropy is given by a simple sum 
of Bekenstein-Hawking entropy Sg and equilibrium entropy of those matter fields. This simple sum of equilibrium 
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entropies is obtained in a general relativistic framework using the Euclidean path-integral method for a black hole 
spacetime and equilibrium matter fields on it. Extending those equilibrium results, we assume that the total entropy 
Stot of the whole system which consists of evaporating black hole and nonequilibrium matter fields is also given by 
a simple sum, 

Stot Sg + Sjn , (5-1) 

where Sg is the black hole entropy given in Ea. (j2.ip and Sm is the nonequilibrium entropy of the general self- 
interacting matter fields including Hawking radiation. Here we consider a quasi-equilibrium regime of black hole 
evaporation process and use the equilibrium Bckcnstein-Hawking entropy for black hole entropy Sg. The GSL 
requires the following inequality during the black hole evaporation process, 

dStot>0. (5.2) 

More discussions on the validity of additivity (|5.ip and of the existence of a well-defined nonequilibrium entropy 
Sm for an arbitrary matter are given later in Subsec l5.3.2l 

The problem is how to deal with the nonequilibrium entropy 5*^ of matter fields. However, because a general 
definition of nonequilibrium entropy is not formulated, the existing proofs of GSL consider an equilibrium between 
a black hole and a heat bath surrounding it, or assume the existence of a well-defined nonequilibrium entropy of 
arbitrary matters [10l[22]. The equilibrium settings succeeded in proving the GSL, but unfortunately the physical 
essence of GSL remains veiled. To understand the veil over GSL, it is important to know explicitly the situation 
considered in the equilibrium settings; the existing proofs of GSL consider self- interacting matter fields as Hawking 
radiation and any other matters around black hole. Then, shift from equilibrium to black hole evaporation, we can 
recognize there are three physical origins of GSL: 

Origin of GSL (a): Self- interactions of the matter fields of Hawking radiation and around black hole. These 
interactions are collision of composite particles and self-gravitation of matter fields. The self-interactions 
causes a self-relaxation of matter fields and produce their entropy. Self-interacting matter fields have a 
positive entropy self-production rate. 

Origin of GSL (b): Gravitational interaction between black hole and matter fields. This interaction consist of 
curvature scattering, gravitational redshift and so on. The gravitational field around black hole works as if 
a virtual medium on which matter fields propagate, then the composite particles of matter fields interact 
with the virtual medium to result in a relaxation towards some equilibrium state. Therefore the gravitational 
interaction between black hole and matter fields produces matter entropy as well as the origin (a). A positive 
entropy production rate of matters results. 

Origin of GSL (c): Increase of black hole temperature Tg along black hole evaporation due to the negative heat 
capacity Cg < given in Eq. (|2.2p . becomes one of origins of GSL as follows: When dEg < due to the 
evaporation, Tg increases since dTg = dEg/Cg > 0. Then, because the matter fields of Hawking radiation are 
ordinary matters of positive specific heat, we find that, the more evaporation process proceeds, the more matter 
entropy the evaporating black hole radiates out in Hawking radiation. This is not the entropy production 
inside the matter fields during propagating outside the evaporating black hole like origins (a) and (b), but the 
growth of the entropy emission rate by black hole along its evaporation. Because the negative heat capacity 
is a peculiar property of self-gravitating systems, this origin (c) is a self-gravitational effect of black hole on 
its own thermodynamic state. 

In the existing proofs of GSL, all of these origins are included and it remains unclear which of these dominates over 
the others. If the GSL would be proven by considering a situation keeping one of them and discarding the others, 
then we can conclude that the one kept is the physical essence which guarantees the validity of GSL. This section 
aims to reveal the origin (c) is the physical essence of GSL. 

To do so, we utilize the NE model. To pick up the origin (c), we remove the heat bath from NE model, and 
put the black hole in an empty and infinitely large flat spacetime. The black hole is bared in this situation. We 
consider this situation throughout the present section, and call it the hare NE model. 

The bare NE model includes the self-gravitational effect of black hole on its own thermodynamic state through 
equations of states (|2.ip . but ignores the self-interactions (including self-gravitational interaction) of radiation fields 
and the gravitational interaction between black hole and radiation fields. This means that the bare NE model 
includes only the origin (c) and ignores the so-called grey body factor. 

Here it should be pointed out that, while the quasi-equilibrium assumption is also valid for the bare NE model 
for classical size black hole evaporation, but the fast propagation assumption breaks down since the radiation fields 
spread out into an infinitely large space. Hence, we have to take the retarded effect on radiation fields into account. 
However because of the quasi-equilibrium assumption, we ignore the special relativistic Doppler effects due to the 
shrinkage of black hole surface. 
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Here one may show an example as an objection: When an inter-stellar gas collapses to form a star, the self- 
gravitational effect of that gas decreases its entropy. Then the origin (c), self-gravitational effects of black hole, 
may not be the essence of GSL. However as discussed later in Subsec IF.3.31 this objection is not true of the black 
hole evaporation. 

Because the bare NE model does not include the origins (a) and (b), the radiation fields have zero entropy 
production rate during propagating in an empty and infinitely large space. Therefore we find 

dStot > dSNE , (5.3) 

where Stot is the total entropy with general self-interacting matter fields given in Ea. (|5.ip . and Sne is the total 
entropy of the bare NE model. Under the assumption of simple sum for nonequilibrium entropies considered in 
Eq. ljS.ll) . Sne is decomposed as 

Sne ■— Sg + Srad , (5.4) 

where Srad is the nonequilibrium entropy of radiation fields propagating in an empty and infinitely large flat 
spacetime. Therefore, if an inequality dSNE > holds, then the GSL dStot > follows and we can conclude the 
origin (c) is the physical essence of GSL. Next subsection shows dSNE > holds. 



5.2 Time evolution of total entropy Sne 

In this subsection, we calculate a time evolution of SNEit) to show the inequality dSNE > along black hole 
evaporation. Here the time t corresponds to a proper time of a rest observer at asymptotically flat region if we can 
extend the bare NE model to a full general relativistic model. In the foUowings, we obtain explicit forms of Sg{t) 
and Srad{t) as functions of t, then calculate the time evolution of SNE{t) under the assumption of simple sum given 
in Eq.fOll. 



5.2.1 Sg as a function of time 

Time evolution of black hole radius Rg is given by the Stefan-Boltzmann law, 

^^-a'T^Ag, (5.5) 

where a' is the generalized Stefan-Boltzmann constant and Ag := 4 7r_Rg is the surface area of black hole. This 
together with Eq. (|2.ip gives 

where Rq := Rg{0) is the initial radius, and it is assumed that the emission of Hawking radiation starts at t = 0. 
This Rg{t) leads the time evolution of thermodynamic states of black hole and radiation fields. Ea. (|5.6p gives the 
evaporation time (life time) of black hole in the framework of the bare NE model, 

i?3 

tempty 128071"— . (5.7) 

Using Eqs. (|5.6p and (|2.ip . we obtain Sg{t) as a function of time. 

Here let us consider about the quasi-equilibrium assumption. To validate this assumption, a sufficiently slow 
evaporation is required. Then we discuss t^pread and w, where tspread is a time scale of radiation fields to spread 
out into infinitely large space, and f is a shrinkage speed of black hole radius. For the first consider the time scale. 
tspread is givcu by a particle of radiation fields traveling across the size of black hole, 

tspread '■— Ro ■ (5-8) 

This tspread 

gives a typical time scale of radiation fields to spread out into the empty and infinitely large space. 
If t empty is longcr than tspread, wc Can consider the black hole evaporation proceeds slowly. Hence to validate the 
quasi-equilibrium assumption, we consider the case satisfying the following inequality, 

A := = 12807r:^ > 1 . (5.9) 

tspread ^ 

For classical size initial condition Rq > 1 together with Ea. p.3p . this requirement (|5.9p is relevant. For the second 
consider the shrinkage speed of black hole radius. 



V := 



dRgit) 



dt 



R 



^ ^empty ^gij')'^ 



(5.10) 
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where Eq. (|5.6|) is used. To validate the quasi-equilibrium assumption, v should be slow enough, v < 1. This gives 



This requirement together with Eq. (ll.3|) means that, because of N/lBn/AS < 1, the quasi-equilibrium assumption 
is valid at least until one Planck time before evaporation time tempty 



5.2.2 Srad as a function of time 

We will obtain the nonequilibrium entropy S^ad as a function of t by applying the steady state thermodynamics 
for radiation fields. But before proceeding to that calculation, the retarded effect on radiation fields has to be 
introduced into the steady state thermodynamics, since the fast propagation assumption breaks down in the bare 
NE model. In this subsection, firstly we construct a modified distribution function for composite particles of 
radiation fields, then obtain Srad by substituting that distribution function into formulae (|3.3[) and (|3.7p . Because 
of the quasi-equilibrium assumption, we ignore the special relativistic Doppler effects due to the shrinkage of black 
hole surface. 

Hereafter r denotes the areal radius from the center of black hole. To take the retarded effect into account, 
consider a particle of radiation fields emitted by black hole at time t and reaches a spatial point of r at time t (> t). 
The emission time t depends on not only the coordinates (t, r) but also an angle 9 between radial direction of r and 
momentum of the particle under consideration (see left panel in Fig llOp . This emission time t(t, r, 9) is obtained as 
a root of the equation, 

Rg{if^{t-iy+r^-2{t-i)rcos9. (5.12) 

This is the equation of degree six about t, and an appropriate root as the emission time is the maximum root in 
range, < t < t. Although another root may exist in this range, however a non-maximum root corresponds to a 
particle emitted at point q in left panel in FigdO] which is obviously unphysical. The other four roots of Eq. ()5.12p 
may be of complex valued. 

Furthermore the angle 9 has an upper bound, < 9 < 9mit,r). To find an explicit expression for 9m, look 
at a particle of radiation fields emitted by black hole at the initial time t = 0, which we call an initial particle. 
Obviously, we have 9m — for t < r — i?o, since no particle has been reached a spatial point of radial distance 
r. The initial particles emitted in radial direction form the boundary of region filled with radiation fields. The 
initial particles emitted in off-radial directions propagate behind the initial particles emitted in radial directions. 
Therefore, the initial particles reach a point of radial distance r > i?o in a time interval, r — Rq < t < — Rf,. 
Within this time interval, the upper bound 9m is given by initial particles, which is obtained with setting i = in 
left panel in FigfTUl 

cos6l„,(i,r) = ^— , for r > i?o and r - i?o < t < -y/r^ - i?2 . (5.13) 

For t > ^7-2 - Rl, the upper bound 9m is not given by any initial particle but by a particle emitted at point b at 
time tm shown in right panel in Fig llOl Then we find 

cos 9m{t, r) , for r > Rq and ^Jr^-Rl<t, (5.14) 

where the time tmit,r) is a real valued root of the equation, 

= Rg{tm)^ + it-tmf ■ (5.15) 

This is the equation of degree six about tm- The appropriate root for tm should be in the range, < tm < t. This 
root may be degenerated, since the trajectory of this particle is tangent to the sphere of radius Rg(tm), and the 
other four roots may be of complex valued. Finally turn to a point of radial distance r < Rq. It is obvious for this 
point that the upper bound 9m is also given by formula (|5.14p . 

cos 9m{t,r) = , forr<i?o. (5.16) 

r 

From the above, the distribution function of radiation fields is given as 

for 9 < 9m{t,r) 

(5.17) 

, for9>9mit,r) 



d{t,r;uj,9) 



exp 



u/f 
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Figure 10: Left panel shows retarded effect on radiation fields. Right one shows the upper bound of 9 ior r < Ro, 
or r < Rq and t > -^/r^ — i?^ . 



where T = Tg{t{t,r,9) ), w = and the signatures "— " and "+" are respectively for bosons and fermions. Here 
we note, while (i, r)-dependence in d{t,r;uj,0) expresses a spacetime dependence, (w, 0)-dependence expresses a 
dependence on momentum p of particles of radiation fields. 

The nonequilibrium entropy of radiation fields Srad is obtained by substituting Eq. (j5.17p into Eas. (j3.3p and (j3.7p . 



t+B.o 

Srad{t) = I dr4:Trr'^Srad{t,r) 



where 



Srad{t,r) : = 



N 



dy ■ 



(5.18) 



(5.19) 
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where the integrals in Eg. p. lip are used, and we set y := cos0, ym ■= cos Om and gp;i = rif, and Uf for bosons 
and fermions respectively. Since the emission of radiation fields starts at t = 0, the radiation fields fill the space in 
range, 

Rg{t) <r <t + Ro. (5.20) 



5.2.3 Total entropy 

Now we obtain the total entropy Spf e of the bare NE model under the assumption of simple sum as in Eq. ([5 

SNE{t) '■= Sg{t) + Srad{t) , 



(5.21) 



where Sg{t) — tt Rg{t)^ is given by Eq. (l5.6p . Srad{t) is by Eq. (|5.18p and time t corresponds to a proper time of a 
rest observer at asymptotically flat region if we can extend the bare NE model to a full general relativistic model. 
Analytic proof for dSNE/dt > is difficult. So we try to show it numerically. To do so, normalize Sne as follows. 



r := 



r 



T^NEir) 



SNE{t) 



Sne{0) 

Then the normalized total entropy is 



, Crrad{T,x) 



Sradjt, r) 

Sne{0) 



Ml 

Rr, 



(1-r) 



1/3 



-1/1 pX T-\-l pi 

Ene{t) = Xgirf + — dx 



XgifiT,x,y)r' 



(5.22) 
(5.23) 

(5.24) 



where A is given in Eq. (|5.9p . Because Satb does not explicitly depend on Rq (or N), then, once the value of A is 
fixed, the choice of Rq (or N) is arbitrary with adjusting the value of N (or Rq) appropriately to match with A. 
We assume A > 1 as mentioned in Eq. (|5.9p . 

Numerical results are shown in figs[TT] and [121 Parameters Rq and N are chosen so that inequality (|5.9p holds 
and time interval < r < 0.999 is included in range (|5.1ip . Those numerical calculations are carried out in that 
interval < t < 0.999 using Mathematica version 5.2. 

Figln] shows time evolution of total entropy Y^ne for < r < 0.999 and A = 10, 50, 100, 10^, 10*, 10^ and 10^. 
The plotted curves in Fig[Tl]are converging as A increases, and those for A > lO'^ are almost coincident What we 
can find with this figure is as follows: 



A technical comment; since calculations by Mathematica takes a very long time for double integral, Fig llll plots Ejvb for r = 
0, 0.1, 0.2, •• ■ ,0.9, and 0.999 for each A. However any singular behavior (oscillation, divergence and so on) of SjvB is unexpected. 
Therefore, in order to find monotone increasing nature of Sjvb, it is enough to show plots of some representative points like Fig llll 
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Figure 11: Time evolution of normalized total entropy Y^ne{t) := SNE{t)/ Sne{0) for < r < 0.999, where r is a 
time normalized by evaporation time tempty The plotted curves are converging as A increases. 




Figure 12: Spatial and temporal evolution of entropy density arad{T, x) and its spherical shell density iirx'^aradiT, x) 
for A = 10"*^, where arad{T, x) := Srad{t, r)/ S]s[e{0) is a normalized spatial density of radiation fields, x := r/Ro is a 
normalized spatial distance and t := t/tempty is a normalized time. The arrows denote direction of null world lines 
(geodesies). 

• Because the plotted curves are converging, it is concluded that the total entropy T,ne is monotone increasing 
for A > 1, and the GSL is well supported as explained in Eqs. (|5.3p and (|5.4p . 

• Since the normalized black hole entropy Xgir)'^ is obviously monotone decreasing, the entropy of radiation 
fields increases monotonously faster than the decrease of Xg. 

• It is suggested rather well that 'Sne{t) gives almost the same curves in E-r graph for all values of sufficiently 
large A > 10'^, and the final value is approximately given by I]Ar£;(0.999) ~ 1.33. 

Left panel in Fig ll2l shows the normalized spatial entropy density for radiation fields (TradiT, x), and right panel 
shows its normalized spherical shell density Airx'^aradi'T, x). Fig[T2]is made with A — 10'*, however the same behavior 
is obtained for all values of sufficiently large A > 10'^ as far as the author checked. The arrows in FiglT^] denote a 
direction of null world lines of particles of radiation fields. What we can find with this figure is as follows: 

• Tracking the graphs in Fig l 121 along a null word line, we recognize the following: Although the spatial entropy 
density of radiation fields decreases (left panel), the spherical shell density on each spherical shell remains 
constant while it spreads out into an empty and infinitely large flat spacetime (right panel). This is consistent 
with the absence of self-relaxation of radiation fields and gravitational interaction between black hole and 
radiation fields. 

• The entropy of radiation fields at black hole surface (right panel) , which equals the amount of matter entropy 
emitted by black hole at each moment, increases monotonously as the black hole evaporation proceeds. This 
means an accelerated entropy emission by black hole. 

Hence we find from figsfTT] and [T^ that the entropy emission rate by black hole increases faster than the decrease 
of black hole entropy. This indicates the physical essence of GSL is neither the origin (a) nor (b) , but the origin (c) 
which is the self-gravitational effect of black hole causing a negative heat capacity and an increase of black hole 
temperature along its evaporation. 
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5.3 Summary and supplementary discussions of this section 

5.3.1 Summary 

Considering the GSL in the context of black hole evaporation, we recognize there are three physical origins of entropy 
increase; (a) self-interactions of matter fields around black hole, (b) gravitational interaction between the matter and 
black hole, and (c) self-gravitational effect of black hole. The origins (a) and (b) give positive entropy production 
rates inside the matter fields. The origin (c) appears as an increase of black hole temperature due to the negative 
heat capacity and gives an increasing entropy emission rate by black hole. Then we consider the bare NE model 
constructed by removing heat bath from the NE model. The bare NE model describes a black hole evaporation 
including only the origin (c) and discarding the others (a) and (b). Applying the steady state thermodynamics to 
radiation fields in the bare NE model, we can calculate explicitly a time evolution of total entropy Sne, and find 
Sne increases monotonously. This denotes the GSL holds as explained in Eqs. (|5.3p and (|5.4p . Hence we conclude 
as follows: 

Physical essence of GSL: The self-gravitational effect of black hole which appears as an increasing entropy 
emission rate by black hole guarantees a validity of GSL. The entropy emission rate by black hole increases 
faster than the decrease of black hole entropy. The entropy production inside the Hawking radiation and the 
other matter fields around black hole is not necessary for a validity of GSL. 

The increase of entropy emission rate is shown by an accelerated increase of entropy density Srad at black hole 
surface shown in Fig [T^ and causes the accelerated increase of total entropy Sne shown in FigfTTl 

Furthermore we find another interesting result from FigUT] An inequality A > 10^ together with Ea. (|1.3p 
corresponds to black holes of initial radius Rq > 5. Then FigfTT] shows that the normalized total entropy I]jv£;(t) 
for black holes of initial radius i?o ^ 5 evolves in almost the same increasing fashion and reaches almost the same 
final value ~ 1.33. On the other hand inequality (|5.3p denotes Sne is the lowest estimate for total entropy of 
a full general relativistic black hole evaporation including self-interaction of matter fields around black hole and 
gravitational interaction between black hole and the matters. Hence we find a result as follows: 

Lower bound for growth of the total entropy: For black hole evaporation in a full general relativistic frame- 
work including self-interaction of matter fields and gravitational interaction between black hole and the mat- 
ters, the final value of total entropy Stot (not its lowest estimate Sne) should be larger than 1.33 x ^^(O). 

Exactly speaking our analysis based on the bare NE model covers only a semi-classical evaporation stage until one 
Planck time before the time t^mpty as discussed in Eq. (|5.1ip . However it is reasonable to expect that the total 
entropy after the end of quantum evaporation stage is greater than the total entropy just at the onset of quantum 
evaporation stage, since a well-defined entropy has to be of non-decreasing. Therefore the above result on the lower 
bound of total entropy should be true of the end state of quantum evaporation stage. 

5.3.2 Supplementary discussion 1 

We discuss about two key assumptions remained to be proven, the existence of a well-defined nonequilibrium entropy 
Sm of an arbitrary self- interacting matter fields and the additivity of equilibrium and nonequilibrium entropies (jS.ip . 
These two assumptions seem to be reasonable as follows: 

On the existence of a well-defined nonequilibrium entropy S'„j, we refer to the present status of study on 
laboratory systems which have self-interactions but not self- and external-gravitational interactions. For example, 
consider a laboratory system is in a nonequilibrium state which is far from an equilibrium but whose heat flux is not 
extremely strong. Then the extended irreversible thermodynamics [4] gives a well-defined nonequilibrium entropy 
fiux up to the second order in expansion by the heat flux. On the other hand, the evaporation time iempty of a 
black hole has a very long time scale as explained in Ea. (|5.9p . The black hole evaporation proceeds so slowly that 
the nonequilibrium state of self-interacting matter fields around black hole is not extremely far from an equilibrium. 
Therefore we can expect that the expansion of their thermodynamic quantities by the heat flux up to the second 
order is a good approximation. Hence it is reasonable to assume the existence of a well-defined nonequilibrium 
entropy of an arbitrary self-interacting matter fields according to extended irreversible thermodynamics. 

Next consider about the additivity (jS.ip . As mentioned in the previous paragraph, the nonequilibrium entropy 
has already been defined well up to the second order in expansion by heat fiux J4j . Furthermore that nonequilibrium 
entropy satisfies the additivity as in ordinary equilibrium thermodynamics. Therefore we may expect that, when 
an equilibrium entropy satisfies the additivity even under the effect of gravity, the nonequilibrium entropy in the 
framework of extended irreversible thermodynamics also satisfies the additivity. Hence, because of the additivity of 
equilibrium entropies of black hole and matter fields j23j , it seems reasonable to assume the additivity (|5.ip . 
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5.3.3 Supplementary discussion 2 



Finally we try to answer an objection: When an inter-stellar gas collapses to form a star, the self-gravitational effect 
of that gas decreases its entropy Then the black hole evaporation under the self-gravitational effect of black hole 
may not result in an increase of total entropy. 

When an interstellar gas collapses to form a star, it is commonly believed that there arises the increase of 
net entropy of total system which consists of the collapsing gas and the radiated matters from the gas. However 
the self-gravitational effect of the collapsing gas causes the decrease of entropy of the collapsing gas. It is briefly 
explained as follows: Since the pressure of that gas at its surface is zero, the loss of energy of collapsing gas par a 
unit time A_E due to energy emission is actually the loss of heat due to the first law of thermodynamics, 

AE^TAS (5.25) 

where T is the temperature of collapsing gas, and AS is the "loss" of entropy of collapsing gas par a unit time. The 
energy loss AE is the minus of the luminosity L of collapsing gas, L = —AE. Then, with the assumption of local 
mechanical and thermal equilibrium of collapsing gas at each moment of its collapse, relation (|5.25p is rewritten 
more exactly to the following form [5], 

L = y^^^-TAS. (5.26) 
where U is the total internal energy of collapsing gas, and Vl is the total self-gravitational potential given by 

n= r dm (~^\ , (5.27) 







r(m) 



where G is Newton's constant, M is mass of collapsing gas, and the radial distance from the center of gas r(m) is 
expressed as a function of mass m inside a sphere of radius r. 

From the above, we recognize that the radius r(m) becomes smaller as the gas collapses, then the self- 
gravitational potential decreases to result in the radiation of energy L > and the loss of entropy AS* < 0. 
This is a similar phenomenon to the evaporation of black hole itself with decreasing its entropy dSg < 0. On 
the other hand, if we consider not a collapsing but an expanding self-gravitating gas, the radius r(m) increases to 
result in the increase of entropy AS > 0. Here we point out that the Hawking radiation in black hole evaporation 
corresponds not to a collapsing gas but to an expanding gas. If we consider the radiation fields of Hawking radiation 
including their self-gravitational effect, the entropy of radiation will increase during spreading out into an infinitely 
large space. This is just the origin (b) of GSL. Hence it seems that the entropy of self-gravitating matter fields 
of Hawking radiation is larger than that of a non-self-gravitating Hawking field. This implies Sm > Srad, and the 
objection mentioned above is not true of the black hole evaporation process. 

Finally let us recall a statement given in the second paragraph; it is commonly "believed" that there arises the 
increase of net entropy of total system which consists of the collapsing gas and the radiated matters from the gas. 
One of the reasons why it is not proven but "believed" is that there has not been nonequilibrium thermodynamics to 
treat the net entropy. Although we considered the black hole evaporation in this section, a similar method based on 
the steady state thermodynamics will be applicable to a star formation process including radiations from collapsing 
gas. Then "believed" will become "proven" . 



6 Concluding remarks 

Exactly speaking the black hole evaporation is a nonequilibrium process. We used the NE model as a simpli- 
fied thermodynamic model of black hole evaporation, and applied the steady state thermodynamics for radiation 
fields |13j . It is this steady state thermodynamics that enables us to treat the nonequilibrium nature of black hole 
evaporation. In the framework of NE model, we can find a detailed picture of evaporation process with energy 
accretion onto black hole [11], which is summarized in Subsec 14771 Also, by modifying the NE model to describe 
the evaporation process in an empty space ignoring grey body factor for the Hawking radiation, we can find the 
essence of the generalized second law of black hole thermodynamics in the context of black hole evaporation [5T] , 
which is summarized in Subsec lSTSl 

One of the results in Secl4]is that the larger the nonequilibrium region around black hole, the more accelerated 
the black hole evaporation under the influence of energy exchange between black hole and its outside environment 
(heat bath). On the other hand, SecE] considered the bare NE model which is obtained by removing the heat bath 
from the NE model. The bare NE model describes the black hole evaporation in an empty space with ignoring grey 
body factor, which is considered in Subsec l4^ to compare with the NE model. Then one may think that the bare 
NE model is a limiting case of NE model which has the infinitely large nonequilibrium region around black hole, 
and the evaporation time of the bare NE model would be zero. But here recall that the fast propagation assumption 
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breaks down for such limiting situation, and the analysis given in Sec|4]may not be applied to the limiting situation. 
Indeed, as discussed in Subsec H.4.31 the bare NE model seems not describe the limiting situation of NE model. 
The bare NE model is quite different from the NE model, because no energy accretion is in bare NE model while it 
is in NE model. The absence of energy accretion makes the bare NE model a qualitatively different situation from 
the limiting situation of NE model. 

As mentioned SecI^J the NE model is not a full general relativistic model, and ignoring gravitational redshift 
and curvature scattering on radiation fields propagating in hollow region. Towards a general relativistic NE model, 
we have to extend the steady state thermodynamics to its general relativistic version. If we will construct a general 
relativistic steady state thermodynamics for radiation fields, the gravitational effects on radiation fields will be 
included in the NE model. Even when the extended steady state thermodynamics for radiation fields are applied, 
the semi-classical evaporation stage is divided into two stages as shown in SeclH quasi-equilibrium one and highly 
nonequilibrium dynamical one. If we concentrate on the quasi-equilibrium evaporation stage, then the evaporating 
black hole can be expressed approximately by equilibrium solutions of Einstein equation. In that case, we can avoid 
the mathematical and conceptual difficulties for describing evaporating black holes (see the beginning of Sec[T]). 
However, in considering the highly nonequilibrium dynamical evaporation stage in a general relativistic setting, we 
will face those difficulties. 

Next let us give a comment on the lower bound of total entropy of black hole and matter fields, summarized in 
Subsec 15751 Concerning the entropy in any spacetime, the so-called covariant entropy bound conjecture (holographic 
conjecture) is very interesting, which, inspired by some quantum gravitational evidences, conjectures a universal 
upper bound on the entropy of any gravitating system [24| . As far as the present author considers, a weak point 
of the covariant entropy bound conjecture at present seems a lack of evidences related to some nonequilibrium 
thermodynamics and/or statistical mechanics, since the existing evidences seems based on the (local-)equilibrium 
hypothesis. When the quantum gravity will cover the very dynamical and nonequilibrium phenomena, the covariant 
entropy bound conjecture should also be related to some nonequilibrium entropy. Then, since our analysis in this 
article is based on the steady states of radiation fields, some insight into quantum gravity may be extracted by 
combining our result of the lower bound of total entropy with the covariant entropy bound conjecture. But at 
present, the research related with the conjecture is left as an open issue. 

Finally focus our comments on ones given from the point of view of nonequilibrium physics. All of the re- 
sults in this article are obtained by NE model and based on the steady state thermodynamics for radiation fields. 
Nonequilibrium thermodynamic approach may be a powerful tool for investigating the nonlinear and dynamical 
phenomena. While we considered the radiation fields which is of non-self-interacting, however nonequilibrium ther- 
modynamics for ordinary dissipative systems has already been established in, for example, the extended irreversible 
thermodynamics [4] . It is applicable not to any highly nonequilibrium state but to state whose entropy flux is well 
approximated up to second order in the expansion by the heat flux of a nonequilibrium state under consideration. 
It is interesting to consider self-interacting matter field for the Hawking radiation and energy accretion, and apply 
the extended irreversible thermodynamics to those fields. Then we may find a variety of black hole evaporation 
phenomena. Furthermore apart from black hole evaporation, since accretion disks around black hole seem to con- 
sist of dissipative matters in realistic settings, the extended irreversible thermodynamics may give a unique new 
approach to investigate black hole astrophysics. 

Apart from black hole physics, the steady state thermodynamics for radiation fields may be helpful to understand, 
for example, the free streaming in the universe like cosmic microwave background and/or the radiative energy 
transfer inside a star and among stellar objects. Also an example of possible application to a star formation process 
is explained at the end of Subsec l5.3.3l Keeping future expectation of such applications in mind. Subsec 13721 is 
devoted to exhibit a detail of the steady state thermodynamics for radiation fields. 
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